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ABSTRACT 

We present the angular correlation function measured from photometric samples com- 
prising 1 562 800 luminous red galaxies (LRGs) . Three LRG samples were extracted 
O ■ from the Sloan Digital Sky Survey (SDSS) imaging data, based on colour-cut selections 

' at redshifts, z » 0.35, 0.55 and 0.7 as calibrated by the spectroscopic surveys, SDSS- 

LRG, 2dF-SDSS LRG and QSO (2SLAQ), and the AAOmega LRG survey. The galaxy 
samples cover « 7600 deg 2 of sky, probing a total cosmic volume of ~ 5.5 h~ 3 Gpc 3 . 
The small and intermediate scale correlation functions generally show significant 
C\) ■ deviations from a single power-law fit with a well-detected break at ~ 1 h~ x Mpc, 

consistent with the transition scale between the 1- and 2- halo terms in halo occu- 
pation models. For galaxy separations 1 — 20 /i _1 Mpc and at fixed luminosity, we 
see virtually no evolution of the clustering with redshift and the data is consistent 
with a simple high peaks biasing model where the comoving LRG space density is 
constant with z. At fixed z, the LRG clustering amplitude increases with luminosity 
O^l ' in accordance with the simple high peaks model, with a typical LRG dark matter halo 

mass 10 13 — 10 14 /i _1 M Q . For r < 1 /i _1 Mpc, the evolution is slightly faster and 
the clustering decreases towards high redshift consistent with a virialised clustering 
model. However, assuming the HOD and ACDM halo merger frameworks, ^2 — 3 
perccnt/Gyr of the LRGs is required to merge in order to explain the small scales 
clustering evolution, consistent with previous results. 
KnJ [ At large scales, our result shows good agreement with the SDSS LRG result of 

Eisenstein et al. (2005) but we find an apparent excess clustering signal beyond the 
BAO scale. Angular power spectrum analyses of similar LRG samples also detect a 
similar apparent large-scale clustering excess but more data is required to check for this 
feature in independent galaxy datasets. Certainly, if the ACDM model were correct 
then we would have to conclude that this excess was caused by systematics at the 
level of Ait; 0.001 — 0.0015 in the photometric AAOmcga-LRG sample. 

Key words: galaxies: clustering - luminous red galaxies: general - cosmology: ob- 
servations - large-scale structure of Universe. 
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1 INTRODUCTION 

The galaxy two-point function whether in its correlation 
function or power spectrum form has long been recognised 
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as a powerful statistical tool for studying L arge-Scale Struc- 
ture (LSS) of the Universe (|Peebleslll980l h In an isotropic 
and homogeneous Universe, if the density fluctuation arises 
from a Gaussian random process, the two-point correlation 
function, £(r), and its Fourier transform, P(k), contains a 
complete description of such fluctuations. It has been used 
to measure the clustering strength of galaxies in numerous 
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galaxy surveys (see e.g. Groth fe Peebles!! 19771; Shanks et al.l 



Il989l ; iBaugh fe Efstathioul Il993l ; iRatcliffe et all Il99cf ) and 

the observed £(r) is reasonably well represented by a power- 

8 over a large range of scales, 
where ro is approximately 



law of the form £(r) = (r/ro) 1 
« 100 /i _1 kpc - 10 /i _1 Mpc 
5 ft -1 Mpc. 

More recently, large 
have become av a ilable 

2dFGRS lCoiless et alJ l200lh and these surveys provide 
a perfect opportunity to exploit the two-po int function as 
a tool to constrain cosmological parameters jHawkins et al 



galaxy redshift s urveys 
(SDSS lYork eTZl |2000|. 



20031 : [Cole et~al . 2005; Ei senstein et~aT]|2005l : lTegmark et al 
2009 ; iPercival et al.ll2007f ) which in turn provides an excel- 
lent test for our current understanding of the Universe and 
the processes by which the LSS were formed. 

In the past, when galaxy redshift surveys were less avail- 
able, the angular correlation function, w(6), was heavily 
utilised in the analysis of imaging galaxy samples. The spa- 
tial correlation fun ction, ((r), ca n be related to w(8) via 
Limber's equation (|Limberl Il953h . alternatively w(0 ) can 
be in verted to £(r) using Lucy's iterative technique ( Lucvl 
1974), both approaches providing a means to recover the 3- 
D clustering information numerically. Even today, galaxy 
imaging surveys still tend to cover a bigger area of the 
sky and occupy a larger volume than redshift surveys and 
therefore could offer a route to a more accurate estimation 
of the correlation functio n and power spectrum (see e.g. 
IBaugh fe Efstathiodll993l ). One of the disavantages of using 
w{0) is the dilution of the clustering signal from projection 
and hence any small-scale/sharp feature which might exist 
in the 3-D clustering may not be observable in w(6). 

As mentioned above, the correlation function at small 
to intermediate scales can be approximately described by 
a single power-law which also results in a power-law w{9) 
but with a slope of 1 — 7. However with larger sample 
sizes, recent analyses of galaxy distrib utions start to show 
a deviation from a simple power-law (IZehavi et al.l l2005bl : 
iPhleps et ail 12006): iRoss et all I2OO7I ; iBlake et al.l 120081 . see 
also lShanks et al.l(l983l ~ This poses a challenge for a physi- 
cal explanation and understanding of non-linear evolution of 
structure formation. Several authors attempted to fit such 
correlation function using a descr iption of halo model frame- 
work (e.g. ICoorav fc Sheth1l2002r i invoking a transition be- 
tween 1- and 2- halo terms which occurs at m 1 h" 1 Mpc 
where the feature is observed. This distance scale could po- 
tentially be used as a 'standard ruler' in tracking the ex- 
pansion history of the Universe, provided that its physical 
origin is well understood and the scale can be accurately 
calibrated. 

Another feature in the correlation function predicted 
by the standard ACDM model is the 'Baryon Acoustic Os- 
cillations' (BAO). BAO arise from sound waves that prop- 
agated in the hot plasma of tightly coupled photons and 
baryons in the early Universe. As the Universe expands 
and temperature drops below 3000 K, photons decouple 
from the baryons at the so called 'epoch of recombina- 
tion'. The sound speed drops dramatically and oscillatory 
pattern imprinted on the baryon distribution as well as 
the temperature distribution of the photons which approxi- 
mately 13 billions years after the Big Bang revealed as the 
acoustic oscillations in the temperature anisotropies of the 
CMB. The equivalent but attenuated feature exists in the 



clustering of matter, as baryons fall into dark matter po- 
tential wells after the recombination. In recent years, the 
acoustic peak scale in the LS S has been proposed as a 
potential 'standard rul e r' (e.g. Blake fc Glazebrookl 120031 ; 
iGlazebrook et all 120071 : iMcDonald fc Eisenstein! 120071 1 for 
constraining the Dark Energy equation of state (w = p/pc 2 ) 
and its evolution. 

For the BAO approach to the study of Dark Energy to 
yield a competitive result, a lar ge survey of several million 
galaxies is generally required (IBlake fc Glazebrookl [20031: 
Seo fc Eisenstemll2003l ; IParkinson et al.ll2007l : lAngulo et all 



20081 ). A galaxy spectroscopic redshift survey would require 



a substantial amount of time and resources. An alternative 
route which will enable a quicker and larger area covered 
is through the use photometric redshift (photo-z hereafter) 
at the expense of the ability to probe the radial component 
directly. The photo-z errors are usually worse than spec- 
troscopic redshift errors, but this can be compensated by a 
larger survey and deeper imaging. 

The potential of the distribution of Luminous Red 
Galaxies (LRGs) as a powerful cosmological probe has long 
been appreciated jGladders fc Yea l2000l : lEisenstein et al.l 
120011 ). Their intrinsically high luminosities provide us with 
at least two advantages, one being the ability to observe such 
a population out to a greater distance whilst the other is the 
possibility of detecting the small overdensity of the BAO in 
matter distribution at « 100 ft -1 Mpc owing to their high 
linear biafl In addition, their typically uniform Spectral 
Energy Distributions (SEDs) allow a homogeneous sample 
to be selected over the volume of interest. Moreover, the 
strong 4000 A break in their SEDs make them an ideal can- 
didate for the photometric redshift route or even a colour- 
magnitude cut as demonstrated by the success of the target 
selection algorithm of three LRG spectroscopic follow-ups 
using SDSS imaging. In fact, the first clear detection of the 
BAO in the galaxy distributio n came from the analy sis of 
LRG clustering at low redshift (|Eisenstein et al]|2005l ). 

Here we shall present new measurements of the angular 
correlation functions determined from colour selected LRG 
samples. We shall show that this route provides redshift dis- 
tribution, n(z), widths that are close to the current photo-z 
accuracy, with none of the associated systematic problems. 
Indeed, one of our aims is to assess the efficiency of this 
route to BAO measurement compared to a full 3-D redshift 
correlation function. This possibility arises because the n(z) 
width that we obtain is comparable to the ~ 100 ft -1 Mpc 
scale of the expected acoustic peak. 

A similar clustering analysis meas uring w(6) o f LRG s 
with photo-z's has been carried out bv lBlake et al. I (|2008h . 
Equipped with a higher-redshift LRG selection algorithm 
whose effectiveness has been tested with the new LRG spec- 
troscopic redshift survey, the VST-AAfi ATLAS pilot run 



dRoss et al 



IBlake et al 



120081) . our approach is an improvement over 



( 2008h as it probes an almost four times larger 
cosmic volume and we extend the analysis to large scales to 
search for the BAO peak. 



1 This is the well kn own luminosity depen d ant bias as shown ob - 
servationally by e.g. iNorberg etafl <2002ft ; IZehavi et all j2005bf ) 
and is expected in hierar chical clustering cold dark matter uni- 
verse l lBenson et al.ll200lT ), 
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Table 1 . Summary of the properties of LRG samples used in this 
study. 



Sample 


z 


Number 


Density 
(dog" 2 ) 


Magnitude (AB) 


SDSS 


0.35 


106 699 


f»13 


17.5 sC r pctro < 19.5 


2SLAQ 


0.55 


655 775 


«85 


17.5 < i deV < 19-8 


AAQ 


0.68 


800 346 


wl05 


19.8 < i dcV ^ 20.5 



The layout of this paper is as follows. An overview of 
the galaxy samples used here is given in iJ3] describes the 
techniques for estimating the angular correlation functions 
and their statistical uncertainties. We then present the cor- 
relation results in In iJS] the clustering evolution of these 
LRGs are discussed. We then investigate a possibility of the 
acoustic peak detection in the w{9) from the combined sam- 
ple in SjfJ] Finally, the summary and conclusions of our study 
are presented in SJ7] 

2 DATA 

The galaxy samples used i n this study were selected photo - 
metrically from SDSS DR5 (|Adelman-McCarthv etlufeOQTf ) 
imaging data based on three LRG sp ectroscopic redshift sur- 
veys with z » 0.3 5 , 0.55 and 0.7 (lEisenstein et al.1 l200ll ; 
ICannon et ail l200fj ; (Ross et al.1 l200ct ) . In summary, these 
surveys utilised a crude but effective determination of photo- 
metric redshift as the strong 4000 A feature of a typical LRG 
spectral energy distribut ion (SEP) moves through SDSS 
g, r, i , and z bandpasses l|Fukugita et al. I ll996l : ISmith" et al.l 
I2002D . In each survey, a two-colour system (either g — r ver- 
sus r — iair — i versus i — z) suitable for the desired redshift 
range was used in conjunction with r or i-band magnitude 
to select luminous intrinsically red galaxies. The method 
employed by these surveys has been proven to be highly ef- 
fective in selecting LRGs in the target redshift range. The 
full selection criteria will not be repeated here but a sum- 
mary of the algorithm s and any additional criteria will be 
highl i ghted below (see lEisenstein et al.l 1200 ll ; ICannon et al.l 
120061 ; iRoss et al.1 |200S| for further details). Redshift distri- 
butions, n(z), of the LRGs from the spectroscopic surveys 
utilised in this work are shown in Fig. [T] The LRG sam- 
ples corresponding to the above n(z) have been carefully 
selected to match our selection criteria hence these n(z) will 
be assumed in determining the 3-D correlation functions, 
£(r), from their projected counterparts, w{6), via the Lim- 
ber (|l953T ) equation. 

All magnitudes and colours are given in SDSS AB sys- 
tem an d are corrected for extinction using th e Galactic dust 
map of lSchlegel. Finkbeiner fc Davis! l|l998t ). In this analy- 
sis, we only used the galaxy samples in the most contiguous 
part of the survey, i.e. the northern Galactic cap (NGC). All 
colours describe d below refer to th e differences in 'model' 
magnitudes (see iLupton et al. I l200ll . for a review on model 
magnitudes) unless otherwise stated. 

Hereafter we shall refer to the photometrically selected 
sample (not to be confused with the spectroscopic sample 
from which they are associated) at average redshift of 0.35, 
0.55 and 0.7 as the 'SDSS LRG', '2SLAQ LRG' and 'AAfi 
LRG', respectively. 

© 2011 RAS, MNRAS OOP. [111261 
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Figure 1. Normalised redshift distributions, n(z), of the three 
LRG spectroscopic surveys used as the basis for selection criteria 
in this study. 



2.1 SDSS LRG 

The sample used here is similar to the target sample of the 
recently completed SDSS-LRG spectroscopic survey which 
contains « 100 000 spectra and cover over 1 h~ 3 Gpc 3 . These 
objects are classified as L RGs on the basi s of their colours 
and magnitudes following lEisenstein et all (2001, E01 here- 
after). The sample is approximately volume- limited up to 
z ~ 0.38 and spans out to z ~ 0.5. The selection i s done us- 
ing (g — r) and (r — i) colours coupled with r-band Ipetrosianl 
(| 19761 ) magnitude system. The algorithm is designed to ex- 
tract LRGs in two different (but slightly overlapped) regions 
of the gri colour space and hence using two selection crite- 
ria (Cut I and Cut II in E01) as naturally suggested by the 
locus of early- type galaxy on this colour plane (see Fig. 0. 
The tracks shown in Fig. [2] were construc ted using a spectral 
evolu tion model of stellar populations (|Bruzual fc Charlotl 
l2003h with output spectra mimicking a typical SED of the 
LRGs. The stellar populations were formed at z ~ 10 and 
then evolve with two different scenarios, namely a) passive 
evolution of an instantaneous star formation (single burst), 
and b) exponentially decayed star formation rate (SFR) with 
e-folding time of 1 Gyr. Solar metallicity and lSalpeterl (|l955f ) 
Initial Mass Function (IMF) were assumed in both evolu- 
tionary models. 

We used the same colour-magnitude selection as that 
described by E01 but with additional restriction on the r- 
band apparent magnitudes of the objects, i.e. 7- pc tro ^ 17.5. 
This is due mainly to two reasons, a) to separate out the 
objects with z < 0.2 because Cut I is too permissive and 
allows under-luminous objects to enter the sample below 
redshift 0.2 as also emphasised by E01, and b) to tighten 
the redshift distribution of our sample while maintaining 
the number of objects and its average redshift (see Fig. [3]). 

The selection criteria mentioned above also has another 
star-galaxy separation algorithm apart fr om the pipeline 
PHOTO classification ijLupton et al]|200ll h This was done 
by setting a lower limit on the differences in r-band point- 
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Figure 2. The colour-colour plot of SDSS LRG cut I and 77 
showing their positions on the gri colour plane compared to the 
predicted colour-colour locus (observer frame) of typical early- 
type galaxies as a function of redshift (see text for more details) . 
Each solid circle denotes the redshift evolution of the colour- 
colour tracks at the interval of 0.1 beginning with z = 0.1 (bottom 
left). 
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Figure 3. The number of objects as a function of redshift from 
SDSS LRG spectroscopic redshift survey also shown is the subset 
of Cut I and 77 with additional magnitude cut, r potro ^ 17.5, 
applied. 



200 000 objects which after applying the additional mag- 
nitude cut mentioned above, becomes 106 699 objects and 
results in the sky surface density of about 13 objects per 
square degree. 

2.2 2SLAQ LRG 

The 2dF-SDSS LRG and Quasar Survey (2SLAQ) is the 
spectroscopic follow-up of intermediate to high redshift 
z > 0.4) LRGs from photom etric data of SDSS DR4 
Adelman-McCarthv et al]|2006l ) using the two-degree Field 
(2dF) instrument on the Anglo- Australian Telescope (AAT). 
This survey is now completed and contains approximately 
13 000 bona fide LRGs with over 90 percent at 0.45 < 
z < 0.8 in two narrow equatorial strips covering 180 square 
degrees. The prim ary sample of the survey (Sample 8, 
ICannon et al1l2006l ;C06 hereafter) was selected using (g — r) 
versus (r — i) colours and 'de Vaucouleurs' i-band magni- 
tude (17.5 < idev < 19.8). The colour selection of Sample 8 
is similar to that of Cut II which utilises the upturn of the 
early-type galaxy locus in gri colour plane and hence is im- 
mune against the confusion with the late-type galaxy locus 
at higher redshift (see Fig. 2 in E01) but the scattering up 
in colour of interlopers from lower redshift and contamina- 
tion of M-stars can also affect the accuracy of the selection. 
The latter could be prevented by using a similar method for 
star-galaxy separation as described in the last section but 
in this case we used the i-band magnitude rather than the 
r-band. Following C06, two criteria were used, 



jpsf ~ imodel > 0.2(21 — idev) 

and 



radius, 



dcV(i) 



> 0.2, 



(1) 



(2) 



where radius, 



doV(i) 



is de Vaucouleurs radius fit of the i-band 



photometry. As reported by C06, approximately 5 percent 
of the cool dwarf M-stars is still present in their sample and 
we shall assume this value when correcting for the dilution 
of the correlation signal due to the uncorrelated nature of 
foreground stars and the LRGs. In this work, we only use 
Sample 8 as this provides us with a narrower n(z) and higher 
average redshift than the whole 2SLAQ sample. 

A sample of 655 775 photometrically selected LRG can- 
didates (~ 5 percent stellar contamination) is returned by 
the SDSS DR5 'Best Imaging' database when the Sam- 
ple 8 selection criteria is used in the SQL query from ta- 
ble GALAXY. Objects with BRIGHT or SATURATED or 
BLENDED but not DEBLENDED flags are not included in 
our sample. 



spread function (PSF) magnitudes and model magnitudes 
as most galaxies populate the upper part of tpsf — r mo dci 
space compare to the foreground star of similar apparent 
magnitude. The algorithm has been proven to be quite ef- 
fective (less than 1 percent stellar contamination) for this 
range of redshift and magnitude although Cut II needs a 
more restrictive threshold, rpsF — r mo d c i > 0.5 as compared 
to 0.3 for Cut I. 

In practice, the LRG sample described here can be 
extracted from the SDSS DR5 imaging database using 
the SQL query by setting the flag PRIMTARGET to 
GALAXY_RED. This yields a catalogue of approximately 



2.3 AAD. LRG 

The AAS7-AAT LRG Pilot obser ving run was carried out in 
March 2006 byH oss et all (|2008l , and references therein) as 
a 'Proof of Concept' for a large spectroscopic redshift survey, 
VST-AAf?. ATLAS, using the new AAOmega instrument on 
the AAT. The survey was designed to target photometri- 
cally selected LRGs out to z w 1.0 with the average redshift 
of 0.7. The target sample was ob served in three 2-d egree 
fields including the CO SMOS field (IScoville et alj|2007l ). the 
COM BO-17 Sll field dWolf et alj|200lh . and 2SLAQ d05 
field jCannon et al.ll2006l ). 
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We follow the survey main selection criteria, 19.8 < 
idev ^ 20.5 togeth er with the riz colour cuts as described by 
iRoss et all &2QQ& . In summary, the cut utilises the upturn 
of the early-type galaxy colour-colour locus similar to that 
used by 2SLAQ and SDSS LRG surveys. The turning point 
of the track on the riz colour plane occurs at z = 0.6 — 0.7 
as the 4000 A feature moves from the SDSS r to i band 
whilst this happens at z m 0.4 in the gri case. The selection 
technique has been proven to work reasonably well by the 
observed redshift distribution. This is further confirmed by 
the ongoing AAT-AAS7 LRG project, the down-sized ver- 
sion of the VST-AAfi ATLAS survey, designed to observed 
several thousands of LRG redshifts for photo-z calibration 
and a clustering evolution study. The n{z) (Fig. [TJ used 
in inferring the 3-D clustering information also includes « 
2000 AAf2 LRG redshifts taken during the run in June 2008 
(Sawangwit et al. 2011, in prep.). 

As emphasised by IRoss et aL I (|200Sl ). the stellar con- 
tamination in the sample can be readily reduced to ~ 16 
percent by imposing star-galaxy separation in the z-band 
without any significant loss of genuine galaxies. Although 
the level of contamination could be further reduced by using 
near-infrared photometry, we do not attempt it here as there 
is no infrared survey that covers the entire SDSS DR5 NGC 
sky with similar depth. Therefore we shall use the quoted 
contamination fraction when correcting the measured w(0) 
for the same reason mentioned in H2.2I Since no expr ession 
for star-galaxy separation is given in Ross et all |2008l ). here 
such a procedure is performed using an equation defining the 
dashed line in their Fig. 3, 



Zpsf - Zmodd > 0.53 + 0.53(19.0 — 2 mo dei) 



(3) 



Applying the above selection rules on the 'Best Imag- 
ing' data of the SDSS DR5 yields a photometric sample of 
800 346 high-redshift LRG candidates with the sky surface 
density of approximately 110 objects per square degree. As 
with the 2SLAQ LRG sample, objects with BRIGHT or 
SATURATED or BLENDED but not DEBLENDED flags 
are discarded from our sample. 

3 ESTIMATING w(9) AND ITS ERROR 

3.1 Optimal estimator and techniques 

The two-point correlation function, f(r), measures the ex- 
cess probability of finding a pair of objects separated by 
distance r relative to that expected from a randomly dis- 
tributed process. The joint probability of finding two objects 
of interest (in this case the LRGs) in the volume elements 
SVi and 8V2 separated by a distance r is given by 



SP(r) = n [1 + £(r)] (SVWa 



(4) 



where n is the number space density of the sample. In prac- 
tice, redshift of individual object is required to estimate the 
separation between a given pair. However if such redshift in- 
formation is not available as in this study, the sky projected 
version, w(9), can be used to analyse the clustering property 
of the sample instead. The 2D equivalent of Eq. [4] is 



sp(8) = k [1 + w(o)] sn L sn 2 



(5) 



where K is the surface density of the objects and SP(6) is 
now the joint probability of finding two objects in solid angle 
Sili and <5f22 separated by angle 6. 



Two possible routes for estimating w(9) are the pixeli- 
sation of galaxy number overdensity, 8 g = S n /n and pair 
counting. The pixelisation approach usually requires less 
computation time but its smallest scale probed is limited by 
the pixel size. We choose to follow the latter. To calculate 
w(6) using the pair counting method, one usually generates 
a random catalogue whose angular selection function is de- 
scribed by the survey. The number of random points are 
generally required to be 10 times the number of objects or 
more. This is necessary to reduce the shot noise. Our ran- 
dom catalogue for each sample has ~ 20 times the number 
of LRGs in SDSS and 10 times for 2SLAQ and AAft-pilot 
(see next section for details on how this was achieved). 

We compute w(8) u sing t he minimum variance estima- 
tor of lLandv fc Szafavl (19931) . It is also an unbiased esti- 
mator ( Martinez fc Saar]|2002l ) for the 2PCF as it can be 
reduced to the exact theoretical definition of 2PCF, i.e a 
variance of density fluctuation in Gaussian field, £(r) = 
(<5(x)<5(x + r)). The form of this estimator is 



DD{9) 
RR{6) 



- 2 



(Nrd 

V N 



DR{6) 
RR(6) 



(6) 



where DD(6) is the number of LRG-LRG pairs with angu- 
lar separation within the angular bin centres at 9. DR(8) 
and RR(8) are the numbers of LRG-random and random- 
random pairs, respectively. The N r d/N ratio is required for 
normalisation. N r d is the total number of random points and 
N is the total number of LRGs. We use a logarithmic bin 
width of Alog(6»/arcmin) = 0.176 for 9 = 0.1' to 50' and a 
linear bin width of 20' at scales larger than 50'. 

The uncertainty in the number density of the sample 
could lead to a bias in the estimation of w(ff) when using 
Landy-Szalay estimator especially at large scales where the 
ampli tude is small and hence we also utilise the Hamilton 
(|1993D estimator, given by 

... DD(9)-RR(9) , 

which requires no normalisation. We used the Hamilton es- 
timator to cross-check our wls for each sample and found 
the difference given by the two estimators to be negligible 
in all three samples. 

For the purpose of determining statistical uncertainty 
in our measurement, three methods of estimating the errors 
are considered. The first method is the simple Poisson error 
given by 



(Tp, 



i(<0 



l + w(0) 



(8) 



For the second method, field-to-field error, we split the sam- 
ple into 24 subfields of approximately equal size. These sub- 
fields are large enough for estimating the correlation func- 
tion up to the scale of interest. This is simply a standard 
deviation of the measurement in each subfield from the best 
estimate and is calculated using 



0FtF(0) 



N 



iV 



DRi{6) 
DR{9) 



wi(8)-w(e)} 



(0) 



where N is the total number of subfields, Wi(0) is a mea- 
surement from the ith subfield and w(9) is measured us- 
ing the whole sample. The deviation of the angular corre- 
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lation function computed in each subfield is weighted by 
DRi (8) j ' DR{8) to account for their relative sizes. 

The third method is the jackkmfe resampling. This is a 
met hod of preference in a number of correlation studies (see 



Scranton et~al"1 |2002| ; IZehavi et all l2005al ; iRoss et all 



[20071 ). The jackknife errors are computed using the devi- 
ation of w(8) measured from the combined 23 subfields out 
of the 24 subfields. The subfields are the same as used for the 
estimation of field-to-field error above. w(ff) is calculated re- 
peatedly, each time leaving out a different subfield and hence 
results in a total of 24 measurements. The jackknife error is 
then 



^(0) = ± DRi ' (e) 

i' = l 



DR(8) 



w i ,{e)-w{e)] 



(10) 



where w t i (8) is now an angular correlation function esti- 
mated using the whole sample except the ith subfield and 
DRii (8) j DR(8) is approximately 23/24 with slight variation 
depending on the size of resampling field. 

The w(8) measured from a restricted area are known to 
suffer from a negative offset called 'integral constraint', ic, 
which tends to force the fluctuation on the scales of the sur- 
vey to zero |Groth fc Peeblesll 19771 ). i.e. w cs t (8) = w(8) — ic. 
The integral constraint can be estimated from the random 
pair counts drawn from th e same angular select ion function 
(£Ll]) as the data (see e.g. lRoche fc Ealesl [l999l ) ; 

ZRR(8)w modcl (8) 

ic= — mw) — ' (u) 

where we assume our fiducial ACDM model (see H4.2|) for 
Wmodci- The ic for the SDSS, 2SLAQ and AAfi-LRG samples 
are 4 x 10 -4 , 1.5 x 10 -4 and 8 x 10~ 5 , respectively. These 
are much smaller than the w(8)'s amplitudes in the angular 
ranges being considered in this paper, as expected given the 
large sky coverage of the SDSS data. 

It is well known that the correlation function bins are 
correlated which could affect the confidence limit on the 
parameter estimation performed under the assumption that 
each data point is independent. Comparison of the estimated 
error using the field-to-field and jackknife techniques to the 
simple Poisson error can give a rough estimate of the devia- 
tion from the independent point assumption. This is plotted 
in Fig. [4] which shows that the assumption is valid on small 
scales where Poisson error is a fair estimate of the statistical 
uncertainty. However the same cannot be said on large scales 
where the data points are correlated and the independent 
point assumption no longer holds. At these scales, such sta- 
tistical uncertainty is likely to be dominated by edge-effects 
and cosmic variance. 

Fig. [4] also shows that the errors estimated using field- 
to-field and jackknife method are in good agreement at all 
angular scales except for 2SLAQ and AAS1 samples where 
the jackknife errors are slightly smaller towards the large 
scales but still agree within 10 percent. The errors quoted 
in later sections are estimated using the jackknife resampling 
method. 

The covariance matrix allows the correlation between 
each bin to be quantified and can be used in the fitting 
procedure to de-correlate the separation bins. We calculate 
the covariance matrix from the jackknife resampling using 
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Figure 4. The ratio of jackknife to Poisson and field-to-field er- 
rors on the measurements of w(8). The open diamonds, triangles 
and solid circles give the error ratios of w(9) estimated from SDSS, 
2SLAQ, and AAfl LRG, respectively. 



where w(6j) is the mean angular correlation function of all 
the jackknife subsamples in the jth bin. Note that the dif- 
ference between w(8j) and w(8) estimated using the whole 
sample is negligible. We then proceed to compute the 'cor- 
relation coefficient', rij, defined by 

(13) 



C« = (N- l){[w(6i) - w(8i)] ■ [w{8 3 ) - wifli)]) 



(12) 



Fig.[5]shows the correlation coefficients for the three samples 
which are strongly correlated at the largest scale considered 
and less at small scales confirming the simple correlation test 
using Poisson errors. Note that for the purpose of model fit- 
tings in the large-scale sections ( 94.21 14.31 and [6} where a 
more stable covariance matrix is required, we increase the 
number of resampling fields to 96 sub-regions with approxi- 
mately equal area. The size of these sub-regions are also big 
enough for the largest scale being considered in this paper. 
The correlation coefficients constructed from these 96 JK 
resampling are shown in Fig. lAljfor the three LRG samples. 

We use the fed-trees code (|Moore et al] l200lh to min- 
imise the computation time required in the pair counting 
procedure. The angular correlation function is estimated us- 
ing the method described above and then correct for stel- 
lar contamination which reduce the amplitude by a factor 
(1 — /) 2 , where / is the contamination fraction for each sam- 
ple given in Sj2] 

3.2 Constructing random catalogues 

In order to calculate the angular correlation function accu- 
rately, a random catalogue is required. This catalogue con- 
sists of randomly distributed points with the total number 
at least 10 times that of the data. Each random point is as- 
signed a position in Right Ascension (RA) and Declination 
(DEC). Since our sample spans a wide range in DEC (see 
Fig. [6] for the SDSS DR5 sky coverage), care must be taken 
to keep the surface number density constant assuming the 
survey completeness is constant and uniform throughout. 
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Figure 5. The correlation coefficients, ry, showing the level of correlation between each angular separation bin for SDSS, 2SLAQ, and 
AAf2 LRG (left to right). Note that for each sample we only show r;j up to the angular separation corresponds to 20 Mpc where 
later we shall attempt to fit power-law forms to the measured ui(0)'s. 

and found no significant difference. We follow their approach 
by excluding the overlap regions and re-calculating the angu- 
lar correlation function of our faintest apparent magnitude 
sample, AAf2-LRG, where the issue is expected to be the 
most severe. We found no significant change compared to 
our best estimate using the whole sample. 

3.3 Inferring 3-D clustering 

The angular correlation function estimated from the same 
population with the same clustering strength will have a 
different amplitude at a given angular scale if they are at 
different depths (redshifts) or have different redshift selec- 
tion functions, 0(z). Therefore in order to accurately com- 
pare the clustering strengths of different samples inferred 
from w(8), one needs to know the sample <j>(z). Even if the 
redshifts of individual galaxies are not available, their 3-D 
clustering information can be recovered if the sample red- 
shift distribution, n(z), is known. The equation that relates 
the spatial coherence length, ro, to the amplitude of w(6) is 
usually referred to as Limber's equation. 

Recently, the accuracy of Limber's equation has been 
called into question. This is due to the assumption made for 
Limber's approximation that the selection function, 4>(z), 
varies much more slowly than £(r) in addition to t he flat- 
sky (s mall angle) approximation. It was shown by ISimonl 
(|2007l ) that such an assumption would lead to w(6) being 
overestimated at large angle where the breakdown scale be- 
comes smaller for narrower (j>(z) (see Fig. [7J. Here, we shall 
use the r elativistic generalisatio n of Limber's equation sug- 
gested bv lPhillipps et all (| 19781 ) b ut without the approx ima- 
tion mentioned above. Following IPhillipps et al.l (1 19781 ) for 
the comoving case, 



Figure 6. An equal area Aitoff projection of a random catalogue 
described in £13,21 The red/grey highlighted regions indicate the 
areas where adjacent stripes are overlapped. Note that the shad- 
ing is purely diagrammatic to show the overlap regions and is 
unrelated to galaxy density. 



Only the random points that satisfy the angular selection 
function of the survey as defined by the mask are selected. 

The mask is constructed from 'BEST' DR5 imaging sky 
coverage giver0 in the survey coordinate (A, rf) and stripe 
number. The sky is drift scanned in a stri p parallel to r\ an d 
two strips are required to fill a stripe (|York et alj l200oh . 
Each stripe is 2.5° wide and their centres are separated by 
2.5°. In addition to the 'BEST' sky coverage mask, we also 
exclude regions in the quality 'holes' and regions defined 
as 'BLEEDING', 'BRIGHT_STAR', 'TRAIL' and 'HOLE' 
in the 'mask' table given by the SDSS database. The final 
mask is applied to both our data and random catalogues. 

Note that further away from the survey equator 
(RA2000 = 185°), the adjacent stripes become overlapped 
which account for almost 20 percent of the sky coverage. 
The 'BEST' imaging database only keep the best photome- 
try of the objects which have been detected more than once 
in the overlap regions. At the faint magnitude limit of our 
sample, this could lead to a higher completeness in the over- 
lap region and introduces bias in the estimate d correlation 
function. This issue has also been addressed bv lBlake et al.l 
|2007i ). They compared the measurement from the sample 
which omits the overlap region against their best estimate 



http : //www. sdss . org/dr5 



w(9) 



f °° dz 1 f(z 1 ) J °° dz 2 f(z 2 )gr) 
Uo X dzf(z)] 2 



(14) 



The source's radial distribution, f(z), is simply given by the 
galaxy selection function, 4>{z), as 

f( z )^ x 2 ( z )^ln c (z)^z) (15) 
az 

where \ is the radial comoving distance, n c (z) is the comov- 
ing number density of the sources and r = r(8, zi,z 2 ) is a 
comoving separation of the galaxy pair. We shall assume a 
spatially flat cosmology (see £|4.2[1 hence 



© 2011 RAS, MNRAS 000,[Tt[26l 



8 U. Sawangwit et at 



1.00 




0.01 



Figure 7. The angular correlation function computed using the 
full (dashed-lines) and approximate (solid-lines) Limber equa- 
tion, derived using a power-law, = (r/ro) ' where ro = 
10 h _1 Mpc and 7 = 1.8 with the SDSS LRG n(z) for the thin 
lines and much narrower n(z) (±0.01 centred at 2 = 0.35) for the 
thick lines. 



r = \/x 2 {zi) + X 2 {z2) - 2x(zi)x(22)cos( 



(16) 



Note that Eq. [14] can also be used to relate a non-power- 
law spatial correlation function to w(6) unlike the con- 
ventional power-law approximation of Limber's equation 
(|Phillipps et al.lll978l ). 

Fig.[7]shows w(8) computed using Eq. [14] (dashed lines) 
compared to the conventional Limber's approximation (solid 
lines) for a power- law £(r) with clustering length 10 h^ 1 Mpc 
and 7 = 1.8. The effect of a much narrower redshift dis- 
tribution (thick lines) is also shown where the break scale 
becomes smaller and the power- law slope of w(8) asymptot- 
ically approa ches that of £(r), agreeing with the finding of 
ISimonl j2007r i. We shall use Eq. 1 141 together with the known 
n(z) to infer the 3-D spatial clustering of the LRGs. 

4 RESULTS 

4.1 Power-law fits 

We first look at the angular correlation function measured 
from the LRG sample at scales less than 1° corresponding to 
approximately 20 Mpc where previous studies suggested 
that the spatial 2PCF can be described by a single power- 
law of the form £(r) = (r/ro)~ 7 (typically 7 = 1.8) and a 
single power-law w(8) with slope 1 — 7 is expected (see Fig. 
[7]). However in this study, we find a deviation from a single 
power-law with a break in the slope at ~ 1 h~ x Mpc in all 
three samples (less significant for the SDSS LRG). The mea- 
surement has a steeper slope at small scales (< 1 h^ 1 Mpc) 
and is slightly flatter on scales up to ~ 20 h' 1 Mpc where 
it begins to drop sharply (see Fig. [5] and Fig. [5]). The inflex- 
ion feature at ~ 1 h~ x Mpc has also been reported in the 
spatial and semi-pr ojected, w v (a), correlation function by 
many authors (e.g. IZehavi et al] | 2005al ; iPhleps et al.l I2OO6I ; 
I Ross et al.l 120071 ; iBlake et al.l|2008h and detections go back 
as far as IShanks et al.l "( 1983 ). We shall return to discuss 
these features in the halo model framework (i )4.3p . 

If we first consider w(8) at scales smaller and larger than 
the break point separately, each can be approximately de- 
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Figure 8. The angular correlation function measured from the 
three LRG samples. The solid lines are the projection of best- 
fit double power-law £(r) with ro and 7 given in Table[2]for each 
sample. The break scales occur at approximately a few arcminutes 
depending on the average redshift of the sample. This corresponds 
to a comoving separation of ss 1 h" 1 Mpc (see Fig. [5J. 

scribed by a power-law with a slope of ~ —1.15 (7 = 2.15), 
and « —0.83 (7 = 1.83), respectively. A more detailed anal- 
ysis is performed by fitting a set of models to the measured 
w(8) using a chi-squared minimisation method with the full 
covariance matrix constructed from the jackknife resampling 
(see £|3 . 1 pi . This allows us to quantify the significance of the 
deviation from the single power-law by comparing its good- 
ness of fit to a double power-law. We proceed by calculating 



= J2 Awie^C^Awi 



(17) 



i ,3 = 1 



where N is the number of angular bins, Aw(8i) is the dif- 
ference between the measured angular correlation function 
and the model for the ith bin, and C" 1 is the inverse of 
covariance matrix. 

The single power-law fit is of the form w(6) = 
(6>/6>o) (1 ~ 7) . We also recover the spatial clustering length, 
ro, and its slope through the fitting via Eq. 1141 For a dou- 
ble power-law, the fitting procedure is performed separately 
at the scales smaller and larger than 9b, corresponding to 
~ 1 h~ l Mpc for all three samples (see Fig. [§}. The largest 
scale considered in the fitting for all cases is w 20 ft -1 Mpc 
where a steeper drop-off of w(8) is observed. 

In Fig. [5] the best-fit power-laws for all three samples 
are shown. The summary of the best-fit parameters are given 
in Table [2] Eq. [14] and [17] are then used to find the spatial 
clustering lengths and slopes that best describe our w(8) 
results. The best-fit clustering slopes from ro-7 analysis us- 
ing Limber's equation are in good agreement with that from 
#0-7 and hence we only report the latter in Table [5] If we 
require continuity in the double power-law £(r) at the break 
scale, such a scale can be constrained by the pair of best-fit 
ro-7's for each sample. From Tabled the double power-law 
break for the SDSS, 2SLAQ and AAQ samples are then at 
2.2, 1.9 and 1.3 h~ x Mpc, respectively (see £)5.2I for further 
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Figure 9. The angular correlation function with the best-fit single (red dashed line) and double (blue solid line) power-law for the SDSS, 
2SLAQ and AAf2 LRGs. Lower panels show the fitting residuals for the single (circles) and double (triangles) power-law. 



Table 2. Parameters for the power-law fits to the angular correlation function derived from three LRG samples. The best-fit parameters 
given are defined such that w(9) = (f/fo) 1 7 an d £( r ) = ( r / r o) 7 - The parameters for the best-fit double power-law are given in two 
rows where the 9 < 6 b result is given in the top row. Also given are the corresponding la error for each parameter. 



Sample 


2 


n g 

(h 3 Mpc~ 3 ) 


0o(') 


Single power-law 

7 ro(h~ 1 Mpc) 




«o(') 


Double power- 

7 


law 

ro(/i -1 Mpc) 




SDSS 


0.35 


1.1 x 10~ 4 


1.69 ±0.03 


2.07 ±0.01 


8.70 ±0.09 


16.2 


1.57 ±0.05 
1.05 ±0.09 


2.19 ±0.03 
1.85 ± 0.04 


7.35 ±0.08 
9.15 ±0.16 


2.2 


2SLAQ 


0.55 


3.2 x 10~ 4 


0.87 ±0.01 


2.01 ± 0.01 


7.50 ±0.04 


57.5 


0.83 ± 0.01 
0.60 ±0.03 


2.16 ±0.01 
1.84 ±0.02 


6.32 ±0.03 
7.78 ±0.05 


3.9 


AAS1 


0.68 


2.7 x 10~ 4 


0.57 ±0.01 


1.96 ±0.01 


7.56 ±0.03 


42.8 


0.56 ±0.01 
0.38 ±0.02 


2.14 ±0.01 
1.81 ± 0.02 


5.96 ±0.03 
7.84 ±0.04 


3.4 



discussion of the possible small-scale evolution of £( r ))- By 
assuming the 1 h~ l Mpc break instead of aforementioned 
values, the w{9) is underestimated by ~ 10 percent for the 
SDSS case (less for the other two samples) which is only 
localised to around 9b- The clustering length (single power- 
law), ro, ranges from 7.5 to 8.7 h~ x Mpc, consistent with 
highly biased luminous galaxies. Single power-law fits to the 
data can be ruled out at high statistical significance. While 
the double power-law give better fits to the data than the 
single power-law, their Xred values indicate that such a model 
is still not a good fit to the data, given the small statisti- 
cal errors. Nevertheless, to first order, the double power-law 
fits provide a good way of quantifying the spatial clustering 
strength of the samples via the use of Limber's equation. 

The best-fit slopes at small scales show a sli ght decrease 
with i ncreasing redshift, similar to that found bv lWake et alj 
(|200&j ). The SDSS LRG sample is more strongly clustered 
than the rest as expected. This is simply because the SDSS 
LRG sample is intrinsically more luminous than the 2SLAQ 
and AAQ LRG samples and is not an indication of evolution. 

The galaxy number density (see Table [2| are calculated 



from the unnormalised n(z), assuming the redshift distribu- 
tion from the spectroscopic surveys as desc ribed in jj2] This 
is ga laxy pair- weighted by n 2 (z) (see e.g. iRoss fc Brunnerl 
120091) 

n ^I Az SSv)- 2 ^/I dzn2 ^ (18) 

where f2 bs is the observed area of the sky, x( z ) is the co- 
moving distance to redshift z and c is the speed of light. 
The samples' pair-weighted average redshifts determined in 
the similar manner as n g are consistent with their median 
redshifts and are given in Table [2] 

To this end, we cut back the faint magnitude limit of 
2SLAQ and AAfi LRG's to i dcV < 19.32 and 20.25, respec- 
tively. These cuts are imposed in order to select the sam- 
ples of galaxies whose comoving number densities are ap- 
proximately matched to that of the SDSS LRG. The K + e 
corrected i-band absolute magnitudes of these samples are 
presented in Fig. 1101 We see that their absolute magnitudes 
are also approximately matched. We note that we do not at- 
tempt to match the LRGs' colour of different samples here. 
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Table 3. Properties and the best-fit parameters for double power-law of w(9) measured from the SDSS-density matched samples. 



Sample 


number 


magnitude 


z 


n g 

{h 3 Mpc" 3 ) 


Double power-law 
7 ro(h~ Mpc) 




2SLAQ* 


182 841 


17.5 < i dcV < 19.32 


0.53 


1.2 x 10~ 4 


2.25 ±0.02 
1.80 ±0.02 


6.33 ±0.04 
8.88 ±0.08 


2.1 


A AS!* 


374 198 


19.8 < i dcV < 20.25 


0.67 


1.1 x 10~ 4 


2.20 ±0.02 
1.76 ±0.03 


6.25 ±0.03 
9.08 ±0.06 


1.7 



This would then allow us to roughly constrain the evolution 
of LRG clustering up to z « 0.68 (see [J5J. A summary of 
the properties of these samples and the best-fit parameters 
are given m Tabled The measured w(6)'s are shown in Fig. 

As expected, the amplitudes of the brighter cut 2SLAQ 
and AAfi samples (denoted by 2SLAQ* and AA51* here- 
after) are higher than the original sample. In its raw form, 
w{6) measured from 2SLAQ* increases relative to 2SLAQ 
more than AAS7 relative to AAS1*, due to the narrower red- 
shift distribution of the 2SLAQ* sample. However, if we per- 
form a double power-law fit to these results, the large-scale, 
> 1 h^ 1 Mpc, clustering lengths are very similar and agree 
within w lcr statistical error. To first order these large-scale 
clustering lengths are also consistent with that of the SDSS 
LRG's. We shall investigate the clustering evolution of these 
LRG samples further in ij5] 



4.2 Comparison of the clustering form to the 
standard ACDM model 

We shall compare our w(0) measurements to the predic- 
tions of the standard ACDM model in the linear pertur- 
bation theory of structure growth framework along with 
the non-linear correction. For the theoretical models, we 
first g enerate matter power spectra, usin g the 'CAMB' soft- 
ware (|Lewis. Challinor fc Lasc nbv 2000). In the case of non- 
linear correction, t he software has the 'HALOFIT' routine 
(|Smith et al.ll2003h implemented. Such matter power spec- 
tra, P m (k, z), are then output at the average redshift of each 
sample. The matter correlation function, £ m (r), is then ob- 
tained by Fourier transforming these matter power spectra 
using 



1 

2^2 



2 sin kr 

.(*)* it 



dk 



(19) 



Under the assumption that galaxies trace dark matter 
haloes, the galaxy correlation function, £g( r )> is related to 
the underlying dark matter by the bias factor, b g , via 



h 2 _ Ur) 



6n(r) 



(20) 



Therfore the bias factor is expected to be a function of scale 
unless galaxies cluster in exactly the same manner as the 
dark matter does at all scales. However, at large scales, i.e. 
the linear regime, the bias factor is approxi mately scale- 
independent over alm ost a decade of scales (|Verde et al.l 
12003 : iRoss et al.ll2008h . 

Although we found the clustering lengths and hence the 
amplitude of £ (r) to be very similar for the SDSS, 2SLAQ* 
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Figure 10. Top: The i-band absolute magnitude distri- 
bution of the spectroscopic LRG catalogues. All photom- 
etry is galactic-extinction corre cted using dust map of 
Schl cgel. Finkbeiner fc Davis! dl998t) and K + e corrected to z = 
using the Early-type galaxy templates from iBruzual fc Chariot! 
|2003J). Bottom: The distribution of the absolute magnitude after 
applying a faint limit cut to 2SLAQ and AAS! LRG in order to 
match the comoving number density of the SDSS LRG. 



and AAS1* samples, the evolution in the dark matter clus- 
tering means that the linear bias could be a strong function 
of redshift as we shall see in the next section where we inves- 
tigate the clustering evolution in more detail. The evolution 
of structures in linear theory fr amework is de s cribed by the 
linea r growth factor, D(z), (e.g. [Peebles 19841: [Carroll et all 
Il992f ) such that 

6(r,z) = D(z)6(r,z = 0), (21) 
recall that £(r) = (<5(ri) 5(t2)), where r = |ri — ra|, then 
Z m (r,z) = D 2 (z)Z m (r,0) (22) 
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Figure 11. (a): The angular correlation function measured from the SDSS LRG and the brighter magnitude limit samples drawn from 
2SLAQ and AAf2 sample (symbols). The solid lines are the projection of the best-fit double power-law §(r) with the parameters shown 
in Table [3] For comparison, the dot-dashed and dashed lines are w(9) measured from the whole 2SLAQ and AAQ samples, respectively, 
(b): Same as (a) but now scaled to AAS1 depth and taking into account the relative amplitude due to the different n(z) widths (see text 
for more details). 



The linear growth factor is unity at the present epoch, 
by definition, and decreases as a function of redshift. The 
€,m(r, z) therefore decreases as the redshift increases hence 
given that the number-density/luminosity matched samples 
have similar £ s (r) amplitudes suggests that the bias in- 
creases as a function of redshift. 

We proceed by projecting the predicted £m(j") using Eq. 
1141 Our fiducial models assume a ACDM Universe with 

fi A = 0.73, On = 0.27, /baryon = 0.167, CT 8 = 0.8, 

h — 0.7 and n 3 = 0.95. The linear bias factor is then es- 
timated by fitting the matter w(8) to our measurements 
for the comoving separation of ~ 6-60 h" 1 Mpc, using 
the full covariance matrices. The best-fit linear bias (Xred) 
for SDSS, 2SLAQ*, AAft*, 2SLAQ and AAD. samples are 
2.09T0.05 (1.2), 2.20± 0.04 (0.65), 2.33T 0.03 (0.66), 1.98± 
0.03 (0.53) and 2.07 ± 0.02 (1.2), respectively. The mea- 
sured biases are co nsistent with the resu lts from other au- 
thors. For example. iTegmark et al.l (|2006l ) analysed P{k) of 
SDSS LRG and found b(z = 0.35) = 2.25 ± 0.08 for the 
best-fit erg = 0.756 ± 0. 035 and for our fiducial as this 
becomes b = 2.12 ± 0.12. iRoss et al.1 (|2007n found 2SLAQ 
LRG b = 1.66 ±0.35 u sing r edshift-space distortion analysis. 
iPadmanabhan et al.l (|2007h , using C(l) of SDSS+2SLAQ 
photo-z sample, found that b(z = 0.376) = 1.94 ± 0.06 and 
b(z = 0.55) = 1.8 ±0.04 (assumed as = 0.9), for our fiducial 
a s these are b = 2.18T0.07 and b = 2.02T0.05, respectively. 

Fig. lllb shows the full scaling of of the w(8)'s, account- 
ing for their survey differences. First, the w(9) of the SDSS, 
and 2SLAQ* samples scaled in the angular direction accord- 
ing to their average redshifts and relative to the AAQ* sam- 
ple. The amplitudes are then scaled to obtain a fair compar- 
ison for samples with different redshift distributions. This is 
done by taking the relative amplitudes of the projections of a 
power-law £(r) of the same clustering strength but projected 
through different n(z) widths. Since the observed large-scale 
clustering lengths are very similar, f» 9 h~ x Mpc, the scaled 



w(6)'s in these ranges agree reasonably well. The figure also 
shows the best-fit biased non-linear model for the AAfT 
sample. Our w(8) shapes in the ranges 6 < r < 60 h^ 1 Mpc 
can be described very well by the perturbation theory in the 
standard flat ACDM Universe (see the Xred f° r the best-fit 
bias factor given above). However, at smaller scales the the- 
ory underestimates the clustering amplitude, as expected for 
early-type galaxies. As we shall see in i ]4.3l that the reason 
for this may lie in the details of how the LRGs populate 
their dark matter halo hosts. 



4.3 Halo model fits 

We fit a halo model 
iBerlind fc Weinberdl20ol ; ICoorav 



Peacock fc Smith! l200d : 
fc Shethll2002l ) to our an- 



gular correlation function results. One of the key ingredi- 
ents of the halo model is the Halo Occupation Distribution 
(HOD) which tells us how the galaxies populate dark mat- 
ter haloes as a function of halo mass. Recently, the model 
has been used to fit various datasets as a means to physi- 



sight into their evolution (e.g. White et al. 2007; Blake et al. 


2008]; IWake et al. 


20081; iBrown et al. 20081; IRoss & Brunner 


20091; [Zheng et al. 


[20091). 



Here, we u s e a three-parame ter HOD model (e.g. 
ISeo et all 120081 ; IWake et al.1 |200S| ) which distinguishes 
between the central and satellite galaxies in a halo 
jKravtsov et alj|2004h . The mean number of galaxies resid- 
ing in a halo of mass M is 

{N(M)) = (N C (M)) x (1 + (N S (M))), (23) 

where the number of central galaxy is either zero or one with 
the mean given by 

<iV c (M)>=exp(^^). (24) 
We assume that only haloes with a central galaxy are allowed 
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Table 4. Best-fit HOD parameters. 



Sample 


z 




Mi 


a 


n g 




-^sat 




Y 2 , 
''-red 






(10 13 /i- 1 M Q ) 


{10 13 h- 1 MQ) 




(10- 4 /i 3 Mpc" 3 ) 


(10 ls h- 1 M Q ) 


(per cent) 




SDSS 


0.35 


2.5 ±0.2 


29.5 ±2.5 


1.58 ±0.04 


1.3 ± 0.4 


6.4 ±0.5 


8.1 ± 1.8 


2.08 ±0.05 


3.1 


2SLAQ* 


0.53 


2.2 ±0.1 


27.3 ±2.0 


1.49 ±0.03 


1.3 ± 0.3 


4.7 ±0.2 


7.0 ±0.8 


2.21 ±0.04 


7.7 


AAS1* 


0.67 


2.1 ± 0.1 


23.8 ±2.0 


1.76 ±0.04 


1.2 ± 0.2 


4.3 ± 0.2 


5.7 ±0.7 


2.36 ± 0.04 


10.1 


2SLAQ 


0.55 


1.10 ±0.07 


13.6 ±1.1 


1.42 ±0.02 


3.2 ± 0.5 


3.4 ±0.2 


10.0 ± 1.1 


1.97 ±0.03 


14.2 


AAQ. 


0.68 


1.02 ± 0.03 


12.6 ±1.0 


1.50 ±0.03 


3.1 ± 0.4 


3.0 ±0.1 


9.0 ±0.09 


2.08 ±0.03 


13.6 



to host satellite galaxies. In such a halo, the satellite galax- 
ies ar e distributed following an NFW profile (|Navarro et al.l 
around a central galaxy at the centre of the halo. We 
also assume that their numbers follow a Poisson distribution 
(|Kravtsov et alj|2004 ) with a mean 



<*<-»-(£)" 



(25) 



The NFW profile is parametrised by the concentration 
parameter c = r v i r /r s where r v i r is the virial radius and r B 
is the characteristic scale radius. We assume iBullock et al.l 
(|200lh parametrisation of the halo concentration as a func- 
tion of mass and redshift, 



c(M, z) 



9 



(1+*: 



(£)' 



(26) 



where M* is the typical collapsing mass and is determined 
by solving Eq. [42] with <r(M») = 8 C (0). 

The galaxy number density predicted by a given HOD 
is then 



n s = dMn(M) (N(M)) 



(27) 



where n (M) is the halo mass fun ction, here we use the model 
given bv lSheth fe Lemsonl (| 19991 ) . The effective galaxy linear 
bias can be determined from the HOD; 



6iin = — / dMn(M)b(M) (N(M)) 



(28) 



where b{M) is the ha lo bias as a function of mass, for which 
we use the model of ISheth. Mo fc Tormenl |200ll ) plus the 
improved parameters of iTinker et al.l ( 20051 ) (see t|5.1.2l Eq. 
I41|) . The average mass of haloes hosting such a galaxy pop- 
ulation is then 



M off = — / AM n(M)M (N(M)) 



(29) 



And the satellite fraction of the galaxy population is given 
by 



1 



dMn(M) (iVc(M)) (N B (M)) 



(30) 



The galaxy power spectrum/correlation function can 
then be modelled as having a contribution at small scales 
that arises from galaxy pairs in the same dark matter halo 
(1-halo term). On the other hand, the galaxy pairs in two 
separate haloes (2-halo term) dominate at larger scales, 

P(k) = Pa + P 2h (31) 

The 1-halo term can be distinguished into central-satellite, 



Pcs(fc), and satellite-s atellite, P ss (k), contributions (see e.g. 
ISkibba fc Shethll2009l ); 



Pcs(fc) 



and 



dM n(M)2 (N C (M)) (N S (M)) u(k, M), (32) 



P EE (k) = \ / dMn(M) (iVc(M)) (N S (M)) 2 u(k, M) 2 , (33) 
n i J 

where u(k, M) is the Fourier transform of the NFW profile 
and we have simplified the number of satellite-satellite pairs 
{N S (N S - 1)) to (N S (M)) 2 , i.e. Poisson distribution. 

For 2-halo term, we implement the halo exclusion, 
'nt-matched' , and scale-dependent halo bias, b(M,r), of 
ITinker et al] (|2005h : 

P 2h (fc,r) = P m (k) x 4- 



Mi im (r) 



dMn(M)b(M, r) (N(M)) u(k, M) 



(34) 



where Pm(k) is a non- linear matter power spectrum (see 
i]4.2[) . Mn m (r) is the mass limit at separation r due to halo 
exclusion an d n' s is the restric ted galaxy number density 
(Eq. B13 of ITinker et~ai1 120051). T he scale-dependent halo 
bias is given by ( Tinker et al.ll2005l ) 



6 2 (M,r) = b 2 (M) 



[l + l.m m (r)Y 



(35) 



[l + 0.69C m (r)] 2 - 09 ' 

where £ m is the non-linear correlation function (see N4.2p . 

The galaxy correlation function is then the Fourier 
transform of the power spectrum which can be calculated 
separately for 1- and 2-halo terms. For the 2-halo term, we 
need to correct the galaxy pairs from the restricted galaxy 
density to the entire galaxy population. This is done by 



1 + U(r)=[^) [l + &(r)] 



(36) 



where £2h( r ) is the Fourier transform of Eq. 1341 

We then project the predicted galaxy correlation func- 
tion to w(8) using Eq. Q3] for a range of HOD parameters. 
The best-fit model for each of our sample is then determined 
from chi-square minimisation using the full covariance ma- 
trix. Note that we exclude angular bins corresponding to 
scales smaller than 0.1 /i -1 Mpc because any uncertainty in 
the £(r) model at very small scales, r < 0.01 ft -1 Mpc, can 
have a strong effect on w(0) even at these scales due to the 
projection. The best-fit M m i n , Mi and a and the associated 
values for n s , M e g , F sa t and by ln are given in Table U The 
la uncertainties on the best-fit M m i n , Mi and a are deter- 
mined from the parameter space where Ax 2 ^ 1. For n g , 
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Figure 13. The best-fit HOD models for the SDSS, 2SLAQ*, AAS1* samples (left) and 2SLAQ, AAQ samples (right). These are scaled 
to the AAf2*/AAr2 depth similar to that shown in Fig. 1 1 lb . The bottom panels show the ratios between the best-fit HOD models and 
the measured correlation functions. 




,16 



M halo (h"H} 

Figure 12. The mean number of LRGs per halo as a function of 
mass (solid lines) from the best-fit HOD for the SDSS, 2SLAQ*, 
AAQ* samples (top) and 2SLAQ, AAf2 samples (bottom). The 
central and satellite contributions for each sample are shown as 
the dashed and dotted lines. 

M e g, F sa , t and 6n n which depend on the three main HOD 
parameters, this becomes A\ 2 ^ 3.53. Fig. [T^] shows the 
best-fit HOD for each sample, the coloured solid lines are 
the mean number of LRGs per halo with the central and 
satellite contributions shown separately as the dashed and 
dotted lines, respectively. 

As expected, the LRGs populate rather massive dark 



matter haloes with the masses ~ 10 13 — 1O 14 /i _1 M0. At ap- 
proximately the same redshift, the more luminous samples, 
2SLAQ* and AAS7*, are hosted by more massive haloes than 
fainter samples. Most of the LRGs, > 90 percent, are central 
galaxies in their dark matter haloes, the satellite fraction is 
only 10 percent or less with the the increasing trend towards 
low redshift. This can be explained in the framework of halo 
mergers at lower redshift (see £|5.2,2[) . The best-fit linear bias 
factors for all samples are in excellent agreement with the 
values derived in £14.21 Also the galaxy number density from 
the best-fit halo model is consistent with that derived from 
Eq. [18] (see Tables [2] and |3j) . 

Note that, to first order, our best-fit HODs are compat- 
ible with the measurements from other authors although a 
direct comparison with samples selected differently may not 
be simple. For example, our SDSS sample has similar space 
density (although at h igher redshift, z = 0.35 versus 0.3) as 
the sample studied bv lSeo et ail (|200Sl ). Our Mi/M min and 
satellite fraction are in excellent agreement with their model 
11 (their best-fit A^-body evolved HOD). But their a is some- 
what lower which is caused by the higher as = 0.9 value 
ijWake et al.ll2008l ) and the lower average redshift. Their Mi 
and M m i n are also somewhat higher than our best-fit values 
for the same reason as for a. Another example, our best-fit 
Mi, Mmin, bim and F aa t fo r 2SLA Q* sample are in good 
agreement with IWake et all (|2008h z = 0.55 2SLAQ selec- 
tion, although our values are somewhat higher which may 
be due to our lower galaxy number density, implying that 
our sample contains rarer and more biased objects. 

The best-fit models for w(9) are shown in Fig. 1131 com- 
paring to the data. Both the models and data are scaled to 
account for the projection effect (see £I4.2[I and are plotted 
at the depth of AAf2*/AAf2 sample. We immediately see 
that while the fits at the large scales (r > 3 h~ 1 Mpc) are 
good, the fits at the small scales and at r ~ 1 — 2 h~ l Mpc 
are rather poor especially for the higher redshift samples. 
This is evident in the high best-fit reduced chi-square val- 



© 2011 RAS, MNRAS 000,rfl[26l 



14 U. Sawangwit et al. 



ues in Table [4] Given our small error bars, this may indi- 
cate that a more complicated halo model may be needed, 
e.g. five/six paramete rs HOD, an impro ved halo-exclusion 
model (see Fig. 11 of lTinker et~alll2005h . or different halo 
concentration parametrisation. Another important point to 
note is that the HOD formalism assumes a volume-limited 
sample, which we do not have here. This means that our ob- 
served galaxy number density corresponds to a cut-off which 
evolves with redshift rather than a cut-off in halo mass or 
LRG luminosity. Nevertheless, to first order the HOD fits 
generally describe the shape and amplitude of our measured 
w{9) and we believe that the derived 6n n and M e g are rea- 
sonably robust despite the statistically poor fits. 



5 EVOLUTION OF LRG CLUSTERING AND 
DARK MATTER HALO MASSES 



5.1 Intermediate scales 



We study the LRGs clustering and dark mat ter halo mass 
evolut ion b y employing the method s used bv lCroom et al.l 
(|2005l ) and Ida Angela et all jiool) to analyse their QSO 
samples. We then proceed by considering the small-scale 
clustering evolution in the framework of the halo model. 



5.1.1 Clustering evolution 

In this section, we make an attempt to quantify the cluster- 
ing evolution of the LRGs via the use of the w(0)'s measured 
from the number-density (roughly luminosity) matched sam- 
ples as presented in the last section. We shall first compare 
the result at the intermediate scal es, 1 > r > 20 ft -1 Mpc, 
to the simple long-lived model of iFrvl 1 19961 ). The model 
assumes that galaxies are formed at a particular time in 
the past and their clustering evolution is determined by the 
influence of gravitational potential where no galaxies are de- 
stroyed/merged or new population created, hence preserving 
the comoving number density. In such a model the galaxy 
linear bias is given by 



b(z) 



=i+m^ 



D(z) 



(37) 



and as we saw in H4.2I that £ m (r,z) = D 2 (z)£ m (r, 0), the 
clustering evolution is such that 



6(0) + D{z) 



6(0) 



Mr,0) 



(38) 



We shall also compare the data directly to the linear 
theory prediction for dark matter evolution in the ACDM 
model, £(r,z) oc D 2 (z). In addition, we shall also check 
the stable clustering and no-evo lution (comoving) clustering 
models of lPhillipps et all (| 1978T ) . The stable model refers to 
clustering that is virialised and therefore stable in proper 
coordinates. For a £(r) with r measured in comoving co- 
ordinates, the stable model has evolution £(r) oc (1 + z) 7-3 
and the no-evolution model has £(r) independent of redshift. 
At these intermediate scales, the clustering is unlikely to be 
virialised so the stable model is shown mainly as a reference 
point. From Eq. 1381 the no-evolution mod el represen ts the 
high bias limit of the long-lived model of IFrvl (| 19961 ). The 
stable and comoving models are similar to the long-lived 



model in that they both assume that the comoving galaxy 
density remains constant with redshift. 

In order to quantify the clustering amplitude of each 
sample, we shall use the integrated correlation function in a 
20 h -1 Mpc spher e as also utilised by several authors (e.g . 
ICroom et al]|2005l ; iRoss et al.ll2008l ; Ida Angela et alj|200a ). 
The volume normalisation of this quantity is then given by 



(39) 



?2 ° = 20« ' C(r)r dT 



The 20 Mpc radius is chosen to ensure a large 
enough scale for linear theory to be valid and in our case the 
power-law with 7 « 1.8 remains a good approximation up 
to ~ 20 h~ Mpc. Furthermore, the non-linearity at small 
scales does not significantly affect the clustering measure- 
ments, when averaged over this range of scales. 

The integrated correlation function, £20, approach also 
provides another means of measuring the linear bias of the 
sample. For this, we again assume scale-independent bias 
which is a reasonable assumption in the linear regime. The 
bias measured in this way is given by 



b g (z) = 



£20,: 



(40) 



The mass integrated correlation functions are again 
computed assuming our fiducial cosmological model using 
the matter power spectra output from CAMB. The values 
for £20, m used here are 0.153, 0.126 and 0.112 for z = 0.35, 
0.55 and 0.68, respectively. 

The £20,9 is calculated using the best-fit double power- 
law parameters for each sample. The results are plotted 
in Fig. 114b along with the best-fit linear theory evolution 
(long-dashed line), stable clustering (dotted line), long-lived 
(dashed line) and no-evolution models (dot-dot-dashed line). 
The linear bias factors measured using the £20 approach are 
given in Tableland also presented in Fig. 114b . The bias fac- 
tors determined here are in good agreement with the large- 
scale ACDM (gjj and HOD (Sg3) best-fit models. 

To extend the redshift range, we shall compare our re- 
sults to t he clustering of early -type galaxies in 2dFGRS 
studied bv lNorberg et all (|2002l ) that roughly match the ab- 
solute magnitude of our samples after the K + e correction. 
These are the samples with —21.0 > Mt, — 51og 10 h > —22.0 
and —20.5 > Mb . — 5 log 10 h > —21.5, being compared to the 
SDSS/2SLAQ*/AAfi* and 2SLAQ/AAQ data and denoted 
N02E1 and N02E2 in Table respectively. We proceed in 
a similar fashion to the procedure described above and use 
the author's best-fit power-law to estimate the £20,9 's and 
hence the bias values (see Table [SJ). 

Both luminosity bins can be reasonably fitted by 
the long-lived model. The best-fit models for the Mi — 
51og 10 h = -22.7 and -22.4 samples have 6(0) = 1.93 ± 0.02 
and 1.74±0.02 with X 2 = 7.34 (3 d.o.f) and 4.11 (2 d.o.f) re- 
spectively, i.e. 1.5 — 1.9cr deviation. This is interesting given 
the lack of number density evolution s een in the LRG lu- 
minosity function bv lWake et al.l (120061 ). Nevertheless, it is 
intriguing that such a simple model gets so close to fitting 
data over the wide redshift range analysed here. 

The stable model and the linear theory (with constant 
bias) model rise too quickly as the redshift decreases, ex- 
cluded at > 99.99% confidence. However, the comoving 
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Table 5. Summary of the estimated LRG and 2dFGRS early- 
type galaxy bias factor and Mdmh as a function of redshift and 
luminosity. 



Sample 


2 


Mi 


6 


-Mdmh 






-51og 10 h 




(1O 13 7i- 1 M0) 


SDSS 


0.35 


-22.67 


2.02 ± 0.04 


4.1 ±0.3 


2SLAQ* 


0.53 


-22.69 


2.16 ±0.04 


3.3 ±0.2 


AAS1* 


0.67 


-22.60 


2.33 ±0.03 


3.1±0.1 


2SLAQ 


0.55 


-22.40 


1.91 ±0.03 


2.1±0.1 


AAS1 


0.68 


-22.37 


2.04 ±0.02 


1.9 ±0.1 


N02E1 


RJ 0.1 


-22.68 


1.90 ±0.23 


6.2 ±2.2 


N02E2 


RJ 0.1 


-22.40 


1.66 ±0.20 


3.9 ±1.5 



model also gives a good fit to the SDSS/2SLAQ*/AAft* 
data in Fig. 114a . as expected from the lack of evolution 
shown in Fig. Illb . For this model to be exactly correct it 
would suggest that there was an inconsistency in these re- 
sults with the underlying ACDM halo mass function. More 
certainly, we conclude that the evolution of the LRG clus- 
tering seems ve ry slow. This general conclusion agree s with 
previous work (White et all 120071 : I Wake et all 12008( 1 The 
latter author also only found a marginal rejection of the 
long-lived model from the large-scale clustering signal (1.8a) 
compared to 1.9a here. They found a much stronger rejec- 
tion of a 'passive' evolution model from the small-scale LRG 
clustering and we shall return to this issue in i\5.2\ 



5.1.2 LRG dark matter halo masses 

The large-scale galaxy bias is roughly the same as that of 
the dark matter haloes which is a known function of mass 
threshold. Thus by measuring the clustering of the LRGs 
one can infer the typical mass of the haloes they reside 
in. The procedure empl oyed here is similar to t hat used by 
ICroom et"afl (120051 ) and Ida Angela et all (|2008h to estimate 
the dark matter halo masses of QSOs. 

An ellipsoidal collapse model r elating a hal o bias fac- 
tor to its mass was developed by ISheth et alj (|200lh as 
an improvement ov er an earlier spherical collapse model of 
iMo fe White[(ll996l). In this ana lysis, we shall use the expres- 



si on given m ISheth et all |2001i ) and the revised parameters 



of iTinker et all (|2005l ) which were calibrated to give better 



fits to a wide range of as values for variants of ACDM model; 



6(Mdmh, z) 



= 1 



y/a5 c ( 



y/a(au 2 ) + ^/ab(av 2 ) 1 c 

K 2 r 



(av 2 ) c + b(l-c)(l-c/2) 



(41) 



where a = 0.707, b = 0.35 and c = 0.80. v is defined as 
v = S c (z) /ct(Mdmh, z). 5c is the critical densi ty for collapse, 
and i s given by S c = 0.15(127r) 2/3 f7 m (z) 0055 (jNavarro et all 
Il997l ). The rms fluctuation of the density field as a function 
of mass Mdmh at redshift z is ct(Mdmh, z) = ct(Azdmh)-D(z) 
where ct(Mdmh) is given by 



c(Mdmh) 2 = — x 
2tv 



k 2 P(k)w(kr) 2 dk 



(42) 



linear theory 
stable clustering 
long-lived model 
no-evolution 

<M,> - -22.7 




0.4 0.6 
Redshift 



Figure 14. (a): The LRG £20 measurements as a function of 
redshift and luminosity. The data at z f» 0.1 (stars) are taken 
from the correlation functions of early-type galaxise in 2dFGRS 
llNorberg et al,|[2002l) . Open and solid symbols correspond to the 
samples with median absolute magnitude, Mi — 51og 10 h = —22.7 
(SDSS/2SLAQ*/AAQ*) and -22.4 (2SLAQ/AAQ). The best fits 
for various models are also shown (see text for more details) . The 
lower luminosity data have been lowered by 0.2 for clarity, (b): 
The LRG linear biases as a function of redshift and luminosity, 
comparing to the best-fit long-lived model, (c): The typical mass 
of dark matter haloes occupied by the LRGs as estimated from 
the halo bias function. The dot-dashed lines are the best-fit evo- 
lution model of dark matter halo mass via the merger framework 
llLacev fe Colelll993l) . 



P(k) is the linear power spectrum of density perturbation s 
and w(kr) is the window function, given by (|Peeblesl ll980l) 



w(kr) — 3 



sin(fcr) — kr cos(fcr) 



(43) 



(kr) 3 

for a spherical top-hat function. The radius r can be related 
to mass via 

1/3 

(44) 



= ( 3 -Mdmh \ 
\ 47rpo ) 



where po = ^mPcrit is the present mean density of the Uni- 
verse, given by p = 2.78 x W^Q^h 2 M© Mpc" 3 . Here, 
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we use the transfer func tion, T(k), fitting formula given 
by lEisenstein fc Hul (| 1998T ) to construct P(k), assuming our 
fiducial cosmology (see i )4,2l) . 

The estimated dark matter halo masses of the LRG 
samples are given in Tableland plotted in Fig. I14fc . Note 
that the formalism of estimating dark matter halo masses 
from the galaxy biases used here assumes one galaxy per 
halo and can overestimate the t hreshold mass for a given 
value of bias (|Zheng et al.l [2T307T I. This is particularly true 
when we consider the mass estimated from Eq. [41] as the 
threshold mass, minimum mass required for a halo to host 
at least one galaxy and compare the results derived here to 
M m in from the best-fit HOD ( )]4.3[1 . However, if it is used 
as an estimate for the average mass of the host halo then it 
is under-estimated by «40 percent compared to the value 
given by the HOD due to the one galaxy per halo assump- 
tion. 

Next, we attempt to fit the derived dark matter halo 
masses of these LRGs to the halo merger framework in hi- 
erarchical models of galaxy formation . We use the formal- 
ism discussed by lLacev fc Cold (|l993u to predict the me- 
dian Mdmh of the descendants of virialised haloes at z = 1 
for a given halo mass and fit this to our data. In essence, 
the model gives the probability distribution of the haloes 
with mass Mi at time ii evolving into a halo of mass Mi 
at time ti via merging. Fig. 114b shows the best-fit models 
for the Mdmh evolution estimated in this way. These mod- 
els appear to be good fits to both luminosity bins with the 
best-fit Mdmh (2 = 1) = 2.32 ± 0.07 x 10 13 h' 1 M Q and 
1.47 ± 0.05 x 10 13 h,- 1 M G for the L > 3L* and > 2L* 
samples, respectively. 

The most massive haloes hosting these luminous early- 
type galaxies appear to have tripled their masses over the 
past 7 Gyr (i.e. half cosmic time) in stark contrast to the 
little evolution observ ed in the LRG stel l ar masses over th e 
same period (see e.g. IWake et al.l 120061 ; ICool et al.l l2008h . 
This lack of evolution contradicts the predictions in the 
standard hierarchical models of galaxy formation where one 
expects the most massive galaxies to form late via 'dry' 
merging of many less massive galaxies. However, this comes 
with a caveat that the Mdm h at z ~ is an extrapola- 
tion (assuming lLacev fc Cold (|l993l ) halo merging model) 
of the z = 0.35 — 0.7 measurements and the constraint on 
the Mdmh (z = 0.1) is much weaker than the higher redshift 
results. 



0.10 



r (h 4 Mpc) 



1.0 



5.2 Small-scale clustering evolution 

Finally, we discuss the evolution of the correlation function 
at scales corresponding to r < 1 h^ 1 Mpc. We concentrate on 
comparing the number density matched AAf2* and 2SLAQ* 
samples to the SDSS sample. As can be seen in Fig. Illb . 
while at larger scales the w(9) show amplitudes that are re- 
markably independent of redshift, at smaller scales the high 
redshift AAS7* sample appears to have a lower amplitude 
than the lower redshift surveys. Here, we compare the clus- 
tering in non-linear regime to two clustering evolution mod- 
els, namely stable clustering and HOD evolution models. 
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Figure 15. The small-scale w(8) at z = 0.35 evolved from 
the best-fit double power-law of AAQ* (green dashed line) and 
2SLAQ* (red dot-dashed line) samples, assuming stable cluster- 
ing model. The ratios of the evolved ui(0)'s to the best-fit double 
power-law of SDSS sample are shown in the bottom panel. The 
shaded regions signify la uncertainties in the best-fit models. 



5.2.1 Stable clustering model 

The stable model describes the clustering in the virialised 
regi me and hence stable ( unchanged) in proper coordinates 
(e.g. iPhillfpps et al.|[l97a i. Therefore, assuming this model 
one expect the spatial correlation function to evolve as 
£(r) oc (1 + z) 7-3 , where r is measured in comoving co- 
ordinates and 7 is the power-law slope of the correlation 
function. Fig. [15] shows the small-scale, r < 1 h^ 1 Mpc, 
w(ff) of the SDSS sample plus its best-fit double power-law 
model, comparing to the evolved w(8) from the zl — 0.53 
and 0.67 best-fit models. Their ratios to the z = 0.35 best- 
fit model are shown in the bottom panel with the shaded 
regions represent la uncertainties in the best-fit models. 
We see that the evolved zl — 0.67 stable model under- 
predicts the z = 0.35 w(9) somewhat but otherwise is within 
the lcr regions of each other with the probability of accep- 
tance P(< x 2 ) = 0.827. The agreement between the evolved 
zl = 0.53 and the z = 0.35 is better, P(< x 2 ) = 0.999, 
given that the redshift difference is smaller. Note that the 
stable clustering model over-predicts the clustering ampli- 
tude at r > 1 h~ Mpc which is also observed in Fig. 114b . as 
expected. 

The physical picture that is suggested is that the inflex- 
ion in the correlation function may represent the boundary 
between a virialised regime at small scales and a comoving 
or passively evolving bias ed re gime at larger scales. As n oted 
bv lHamilton et al.1 (19911 ) and jpeacock fc Doddsl l|l996h . the 
small scale, non-linear, DM clustering is clearly expected 
from TV-body simulations to follow the evolution of the viri- 
alised clustering model. However, for galaxies in a ACDM 
context, the picture may be more complicated. 

For example, by comparing the 2SLAQ and SDSS LRG 
redshift surveys using the semi-projected correlation func- 
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tion, IWake et all (|2008l ) have suggested that a passively 
evolving model is rejected, weakly from the large scale evo- 
lution but more st rongly from the evolution at small scales. 
IWake et all (|200Sft interpret the clustering evolution using a 
HOD description based on the ACDM halo mass function. 
Their 'passive' model predicts a far faster evolution at small 
scales than is given by our stable clustering (see Fig. I16[) . 
Our stable model is certainly passive in that it is based on 
the idea that the comoving number density of galaxies is 
in dependent o f redsh ift. However, the passive HOD model 
of I Wake et all (|2008l ) requires only 7.5 percent of LRGs to 
merge between z=0.55 and z=0.19 to reconcile the slow LRG 
density and clustering evolution in the ACDM model. We 
shall see in the the next section if this model can also ac- 
commodate our z=0.68 clustering result while maintaining 
such a low merger rate. 



5.2.2 HOD evolution 



In 95.1.11 we foun d using th e large-scale linear bias that the 
long-lived model (|Frvlll996l ) is only marginally rejected at 
1.5-1. 9*7. This is in good agreement with the similar analy- 
sis of Wake et al.l (|2008l ). However, they argued that if the 
small-scale clustering signal was also taken into considera- 
tion, the long-lived model can be ruled out at much higher 
significance (>99.9 percent). 

Recall that our goodness-of-fit (based on the minimum 
X 2 ) for the halo models is rather poor (see table [4|. This 
may be an indication that a more complicated model may 
be needed, e.g. five-parameters HOD and/or a better two 
halo-exclusion prescription etc., given our small error bars. 
Nevertheless, the HOD fit generally describes the shape and 
amplitude of our measured w(8) between 0.1-40 ft -1 Mpc. 
Therefore, at the risk of over-interpreting these HOD fits, we 
make a further test of the long-lived model by evolving the 
best-fit HODs of the higher redshift samples to the SDSS 
LRG average redshift. 

Following the methods described in (|Wake et al.ll200il 
and references therein), the mean number of galaxies hosted 
by haloes of mass M at later time, zo, is related to the mean 
number of galaxies in haloes of mass m, (N(m)), at earlier 
time, z\, via 



(N(M)) 



AmN{m,M) (N(m)) 



dmN(m,M) (N c (m)) [1 + (N a (m))} 



= C(M) + S(M), 



(45) 



w here N(m, M) is the c onditional halo mass function 
of ISheth fc Tormenl (|2002h which is the generalisation of 



lLacev fc Cold l| 19931 ) results, C(M) and S(M) are the num- 
ber of objects which used to be central and satellite galaxies. 

We then model the central galaxy counts in the low- 
redshift haloes assuming that the progenitor counts in these 
haloes is 'sub-Poisson ' <|Sheth fc Lemsonl 1 19991 : 1 eo et al.l 
l200Sl ; IWake et al.ll20Qg| ) such that 



(N C (M)) = 1 



C(M) 
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Figure 16. The predicted SDSS LRG w(0) from passively 
(/no-merge = 1) evolving the best-fit HODs of 2SLAQ* (zl=0.53, 
red dot-dashed line) and AAf!* (zl=0.67, green dashed line) sam- 
ples. The results when central galaxies from high redshift samples 
are allowed to merge (see text for more detail) are also shown, 
blue dotted and magenta long-dashed lines. The bottom panel 
shows the ratios of the evolved w(9)'s to the SDSS best-fit, the 
shaded regions signify the 1<t uncertainties. 



where A^max — int (M/M m i n ). This model is favoured by the 
IWake et~al . (2008) an alysis and is also seen in the numerical 
models of lSeo et all (|2008l ). The mean number of satellite 
galaxies in the low-redshift haloes is then given by 



(JVc(M)) (N.(M)) = S(M)+f no - 



[C(M)-(N C (M))}, (47) 



(46) 



where /no-merge is the fraction of un-merged low-z satellite 
galaxies which were high-z centr al galaxies. This m odel is 
called 'central-central mergers' in I Wake et al.l l|2008h . where 
the more massive high-z central galaxies are more likely to 
merge with one another or the new central galaxy rather 
than satellite-satellite mergers. 

For the long-lived model, we set / n o- mC rgo = 1. The re- 
sults of passively evolving the best-fit HODs from zl = 0.67 
(AAf2*) and zl = 0.53 (2SLAQ*) to zO = 0.35 are shown 
in Fig. [TS] along with the SDSS best-fit model. At large 
scales (r ^ 5 h~ x Mpc), the long-lived model can only be 
marginally rejected at no more than 2a for the AAf2* case 
and is consistent within la in the case of 2SLAQ*. How- 
ever, if we now consider the small-scale, r < 1 h^ 1 Mpc, 
clustering signal we see from the bottom panel of Fig. [16] 
that the long-lived model becomes increasingly inconsistent 
with the best-fit model at z = 0.35. For r ^ 0.5 h~ x Mpc, 
the long-lived model can be rejected at 99.88 and > 99.99 
percent significance using the evolved 2SLAQ* and AAfi* 
HODs, respectively. The much higher clustering signal at 
small scales is caused by far too many satellite galaxies in the 
low-redshift haloes being predicted by the long-lived model. 
This also results in the higher satellite fractions than ob- 
served; both evolved 2SLAQ* and AAfi* give F sat = 18 ± 1 
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percent at z 
best-fit. 



0.35 compared to 8.1 ± 1.8 seen in the SDSS 
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Next, we ass ume the central-central mergers model 
|Wake et alj l2008h and attempt to match the large-scale 
clustering signal of the evolved HOD from high-z to the 
z = 0.35 best-fit model. As argued by IWake et all l|2008h 
and here that this is more likely to happen than the satellite- 
satellite merging case. The /no-merge parameters in Eq. [47] 
required to give the best matches to the large-scale cluster- 
ing amplitude of the SDSS best-fit is 0.2 and 0.1 for the 
2SLAQ* and AAfi* case, respectively. The new w(9)'s de- 
termined from these models are plotted in Fig.[T5]as the blue 
dotted and magenta long-dashed lines. We can see that the 
zl = 0.67 evolved w(0) at small scales is in excellent agree- 
ment with the SDSS best-fit model. The predicted satellite 
fraction, F sat = 7.8 ± 0.9, is also consistent with the SDSS 
best-fit value. For the zl — 0.53 case, the small-scale cluster- 
ing signal is still somewhat stronger that the SDSS best-fit 
model but otherwise are within la confidence regions of each 
other, and the predicted F sa ,t = 10.5 ± 1.3 is also somewhat 
higher than the best-fit value. The galaxy number density 
is reduced due to these central-central merger by «6 and 11 
percent for the zl — 0.53 and 0.67, respectively. However, 
note that this is 2-3 times smaller than the fractional errors 
of our best-fit n g , ~ 20 percent. 

In order to get a handle on the merger rat es which can 
then be c ompared to the pr evious results of IWhite et al.l 
|2007l ) and lWake et~all |2008h . we follow their method of ad- 
justing the galaxy number density. This is because for this 
type of analysis the galaxy samples at different redshifts are 
usually designed to have the same space density. Whereas 
merging means that the space density of the low-z sample 
must be reduced unless there are new galaxies created via 
merging of the fainter objects which fail to be in the high- z 
sample but become bright enough to be in the low-z sample. 
To account for such an effect by physically removing galax- 
ies in a sample i s rath e r difficult to d o in p racti ce as argued 
by IWake et al] (|2008l ). IWhite et all (|2007f l and IWake et all 
( 2008]) adjusted the mass-scale of the low-z HOD fit by sev- 
eral per cent which reduce the space density and increase the 
clustering signal and hence require lower amount of merging 
of the high-z population needed to match the low-z measure- 
ment. Increasing the /no-merge factor in Eg. 1471 results in a 
higher galaxy number density and clustering signal. There- 
fore there is only one unique solution of mass-scaling and 
merging fraction that will simultaneously match the galaxy 
number density and the clustering signal (at large scales) of 
the evolved and best-fit HODs at low-z. 

We increase the mass-scale of z = 0.35 HOD fit by 12 
(7) percent and allow 60 (50) percent of the zl = 0.67 
(0.53) central galaxies to merge in order to get the matched 
large-scale bias of 2.12 (2.10) and n e = 1.12 (1.19) x 
10" 4 /i 3 Mpc -3 . This yields the merger rate between z=0.67 
(0.53) and z=0.35 of « 6.6 (5) percent, i.e. « 2.8 (3.4) 
percent Gyr -1 . The evolved w(9) divided by the model at 
z = 0.35 with increased mass-scaled HOD fit is shown in Fig. 
1171 As noted earlier, the reduction in the galaxy number den- 
sity is small compared to its best-fit fractional error which 
means that our contraints on these merger rates are rather 
weak. However, to first order the merger rates derived here 
appear to be consistent with the value of 2.4 ± 0.7 and 3.4 
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Figure 17. The ratio of the evolved w(0) to the SDSS best-fit 
model with the HOD mass-scale increased by 12 percent. 



per ce nt Gyr 1 found bv lWake et alj (|2008T ) an d White et al.l 
l|2007h . respectively. 

In summary, the combination of the stable clustering 
and passive evolution model is remarkably close to explain- 
ing the clustering evolution of the LRGs at small and large 
scales. These models are much simpler than the HOD frame- 
work which require an understanding of how galaxies popu- 
late dark matter haloes and how they and their host haloes 
merge. The galaxy long-lived model in the context of halo 
framework is significantly incompatible with the small-scale 
clustering data and requires that « 2 — 3 percent/Gyr LRGs 
to merge in order to explain their slow clustering evolution. 
On the contrary, the stable model requires the comoving 
number density to be constant with redshift. This may sug- 
gest that the simple virialised model may only provide a 
phenomenological fit to the small-scale clustering evolution 
in the context of the ACDM model. 



6 SEARCHING FOR THE BAO PEAK 

Next, we inspect the correlation functions at larger scales 
to make a search for the BAO feature. We first present the 
raw correlation functions in Fig. 118a . Note that the inte- 
gral constraints (see £|3.ip are sub-dominant compared to 
w(9)'s amplitudes at these scales. Each correlation function 
shows a feature at large scales, the most significant detection 
comes from the AAQ sample where the clustering signal at 
120' < 9 < 500' is detected (above zero) at more than 4<r 
significance, P(< x 2 ) = 1 X 10~ 6 (with covariance matrix) 
and 3.5cr significance for 200' < 9 < 500'. 

The question is are these features real or simply due to 
systematic error? (see fl6.ll for a series of systematic tests) . 
Here, we perform a classic scaling test to see if any feature is 
reproduced at the different depths of the three LRG samples. 
Given that the samples have intrinsically different ro (see 
Table [2]), we choose simply to scale in the angular direction 
only. The SDSS and 2SLAQ LRG correlation functions are 
scaled in the angular direction to the AASTs depth using the 
average radial comoving distance of each sample. In Fig. 118b . 
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Figure 18. (a): The angular correlation function of the three LRG samples at large scales. The shaded regions are la JK errors. The 
arrow indicates the expected BAO angular separation in each sample, assuming our fiducial cosmology, (b): Same as (a) but now scaled 
in the angular direction to the depth of the AAC LRG sample. 
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Figure 19. The combined angular correlation function of the 
three LRG samples scaled to the AAf2 depth, comparing the re- 
sults when the SDSS standard (diamonds) and uber- (circles) 
calibration are used. Also shown is the average ficld-to-ficld w(9) 
(asterisks) which represents an attempt to filter out any large 
scale gradients in the SDSS data. 



we see that the scaling agreement of the large scale, m 300', 
features is poor. Although SDSS shows a moderately strong 
peak feature, this is not reproduced at the same comoving 
physical scale in the other two datasets. 

Despite this failure of the scaling test, we now attempt 
to increase the signal to noise ratio by combining the mea- 
surements from the three samples using inverse quadrature 
error weighting. Firstly, the SDSS and 2SLAQ w(0)'s are 
scaled in the angular direction to the depth of the AAQ, 
LRGs (radial comoving distance, x ~ 1737 h~ x Mpc as op- 



posed to « 1451 h" 1 Mpc for 2SLAQ and « 970 h' 1 Mpc 
for SDSS) where their amplitudes and errors are then in- 
terpolated to the AAfi's angular bins (i.e. Fig. 118b ). The 
amplitudes of the scaled SDSS and 2SLAQ w(9)'s are then 
normalised to that of the AAfi sample's at 10'. This in- 
volves lowering SDSS and 2SLAQ amplitudes by 25 and 15 
percent, respectively. The resulting correlation function is 
presented in Fig. [19] with the arrow showing the expected 
position of the BAO peak. Note that due to the relatively 
small statistical errors of the AAQ LRG compared to other 
samples, the w(0) result is dominated by the AAQ sample, 
therefore the possible SDSS peak at « 100 h^ 1 Mpc is not 
evident in the combined sample. There also seems to be an 
indication of an excess out to possibly 200 Mpc (see £|6 . 1 1 
for a robustness test of this excess clustering signal) . 

Using the ubercalibration jPadmanabhan et al.l l2008h 
instead of the standard calibration, we find similar results 
at small and intermediate scales but somewhat lower ampli- 
tude at ~ 100 h^ 1 Mpc although the results agree within the 
la error (see Fig. I19[) . This means the correlation functions 
at small and intermediate scales including the parameters 
derived (e.g. power-law fits, linear biases, dark matter halo 
masses) in the earlier parts are not affected by which cal- 
ibration we use. The biggest difference, although less than 
1.5a, is observed at scales larger than 120 h~ l Mpc and up to 
150 h~ x Mpc where the correlation signal is small and hence 
more prone to possible systematics. The weak dependence 
of w(0) at very large scales on the different calibrations may 
be an indication that this apparent extra peak at f« 300' 
could indeed be a systematics effect. We shall return to this 

in sjrni 

We also tested whether the 200 Mpc excess can 
be eliminated by taking the average w(0) from 15 x 20 
deg 2 subfields. The result, after integral constraint correc- 
tion, is shown in Fig. 1191 The 200 h' 1 Mpc excess persists 
even though there is some change at smaller scales. Given 
the model dependence introduced by the integral constraint 
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Figure 20. The angular correlation functions for SDSS, 2SLAQ and AAf2 samples (left to right), measured with varying dust extinction 
limit, astronomical seeing and different star-galaxy separation (top row), compared to our 'Default' results. Also shown is the effect of 
low galactic latitude region exclusion for each sample (bottom row). Note that for the b > 60°, the sample size is reduced by 60 percent. 
In each case, an arrow indicates the expected position of the BAO peak assuming our fiducial cosmology. 



correction (Eq. 1111) . hereafter, we shall use the correlation 
function of the ubercal sample measured using our normal 
method. 



6.1 Testing for systematic effects 

We have performed a series of tests to check our results 
against possible systematic effects. The tests include exclu- 
sions of high dust extinction and 'poor' astronomical seeing 
regions, an improved star-galaxy separation for the AAQ 
sample and effects of possible contamination by clustered 
stars. 

First, we exclude the regions where the i-band extinc- 
tion is greater than 0.1 mag which discards w 20 percent of 
the data. The results are shown in the top row of Fig. 1201 
For 2SLAQ and AAQ samples, the results appear to be lower 
than the main measurements but otherwise remain within 
lcr statistical errors of each other. Although the amplitudes 
at 6 ^ 220' are somewhat lower than the default AAQ re- 
sult, the excess at 9 > 300' still persists. We then investigate 
the effect of excluding the regions with 'poor' astronomical 
seeing, the limit of l".7 is used following the SDSS 'poor' 
seeing definition which discards « 30 percent of the data. 
The results here are in good agreement with the main results 
with the exception of a few angular bins around 320' of the 
2SLAQ sample where they are somewhat (non-significantly) 
lower than the default measurements. 

Next, we attempt to reduce the stellar contamination 
fraction in the AAQ sample. As a reminder, our default 
(optimised) star-galaxy separation algorithm (see S[2]) leaves 



~ 16 per cent stellar contamination in the sample while los- 
ing genuine LRGs only at a sub-percent level. Here, we im- 
pose a more aggressive star-galaxy separation cut which re- 
duces the contamination level to «9 per cent at the expense 
of nearly halving the number of genuine AAQ LRGs. The 
cut is a combination of the fitted 'de Vaucouleurs' radius as 
a function of ZdcV magnitude and the correlation between 
the 'de Vaucouleurs' and fiber magnitudes in z-band. The 
w(8) measurement for this new AAQ sample after correction 
by a factor of 1/(1 — /) 2 , where / = 0.09 is shown in the 
top-right panel of Fig. [20] This is in good agreement with 
the main results. 

W e test our earlier asssumption see also lBlake et al.l 
l200ct ) that the effect of the stellar contamination is sim- 
ply a dilution of <5 g by (1 — /) and hence the amplitude of 
galaxy-galaxy correlation function by (1 — /) 2 , where / is 
the contamination fraction. We add a sample of red stars to 
the 2SLAQ sample at the 16 percent level, similar to what 
we expect in the AAQ sample. The stars are selected from 
SDSS photometric objects which are classified as 'star'. The 
colour selections have been matched to that of the AAfi- 
LRG sample. The sample is then randomly selected to have 
the number of objects at 16 percent of the 2SLAQ sample. 
They therefore follow the stellar distribution with galactic 
latitude. The w(8) result after correction by 1/(1 — /) 2 is 
shown in the top-middle panel of Fig. [20] and is found to 
be in excellent agreement with the main 2SLAQ result. We 
do not see any evidence of a slope change which may arise 
from a possible clustering of the stellar contaminants at large 
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in terms of a systematic that may produce excess cluster- 
ing at large scales. However, such an effect would predict a 
decrease in galaxy density at low galactic l atitudes and this 
is no t seen in our samples (see Fig. 9 of ISawangwit et al.l 
l2010t ). If anything, the opposite effect is seen in our data 
with an increase in density towards lower galactic latitude 
which may be caused by stellar contamination. Here, we 
have tested our w(9) measurements by successively cutting 
out data at low galactic latitudes. Although the 2SLAQ re- 
sults may show some marginal dependence on the galactic 
latitude cut, the AAS7 results seem reasonably unaffected 
(see Fig. |2Uf D). This may be due to the higher stellar con- 
tamination fr action in AAO sa mple which means that the 
effect seen bv lRoss et al.l (|201 lh may not be directly appli- 
cable to the A AO. sample. 

We conclude that the apparent clustering excess at 
~ 300' in the AAQ. sample appears to be reasonably ro- 
bust against most of the systematic tests we performed here. 
However, one might argue that the weakening of the ex- 
cess signal when i ex tinc > 0.1 regions (~20 percent) are ex- 
cluded and the marginal dependence on the galactic latitude 
cuts of the 2SLAQ results may be taken as some evidence 
for systematic effects. On the other hand, the SDSS-AAS1 
cross-correlation test also tends to limit the size of possible 
systematic errors. 



Figure 21. The auto correlation functions for SDSS (top panel) 
and AAC (bottom panel) samples, comparing to the CCF be- 
tween the two samples (blue solid lines). 



scales, at least for the contamination level expected in our 
sample. 

We apply various minimum galactic latitude cuts on 
the data in order to test for any systematic error. Such sys- 
tematics (if they exist) could be due to the gradient caused 
by galactic dust extinction and/or different stellar contam- 
ination fractions which one might expect to be worse in 
the lower galactic latitude regions. Note that in our default 
datasets ~ 95 per cent of the data are at b ^ 30° . The results 
of applying the galactic latitude cuts of b > 40° , 50° and 60° 
are shown in the bottom row of Fig. [201 Note that with the 
b ^ 60° limit, « 60 percent of the data are discarded. The 
2SLAQ results appear to be marginally dependent on the 
galactic latitude limits. In the AAQ. sample the results are 
in good agreement with the main measurement although the 
b > 60° limit appears to be « la lower in some angular bins. 

Finally, we cross-correlate the SDSS and AAQ, samples. 
The redshift distributions of the two samples are well sep- 
arated with only slight overlap (see Fig[T]). Therefore any 
residual signal in their cross-correlation function, CCF, at 
large scales can be used as an evidence for systematic errors. 
The CCF is shown in Fig. 1211 comparing to the auto corre- 
lation functions of the SDSS and AAQ samples in the top 
and bottom panels, respectively. The CCF has much lower 
signal than the ACF at < 120' and is consistent with zero, 
P{< X 2 ) = 0.997, between 120' < 6 < 500' whereas the 
AAQ w(8) signal is detected at more than 4cr significance 
(see above) in the same angula r rang es. 

We note that iRoss et aU (|201ll ) have suggested that 
there is a systematic effect associated with the area effec- 
tively masked by foreground stars which may be important 



6.2 Model comparisons 

6. 2. 1 Standard A CDM model 

First, we compare the measured angular correlation func- 
tion to the perturbation theory prediction in the standard 
ACDM Universe. To compute the theoretical prediction, we 
proceed in the same manner as described in i]4.2l calculat- 
ing w(8) by projecting £(r) which is a Fourier transform 
of a non-linear P(k). However, here we assume t he be st- 
fit cosmological parameters from lEisenstein et al.l (|2005h , a 
flat ACDM model with Q m h 2 = 0.13, Q b h 2 = 0.024, h = 0.7 
and n = 0.98. And unlike in £14.21 the non-linear modelling 
of the BAO peak using only HALOFIT is not adequate. The 
BAO peak in the correlation function can also be broadened 
(and perhaps slightly shifted) by the non-linear gravity sup- 
pression of the higher harmonics in the power spectrum via 
mode couplin g (iMeiksin et al. 1 19991) . To model such an ef- 
fect, we follow lEisenstein et al. ( 2005T ) and smoothly interpo- 
late between the linear power spectrum and the 'no-wiggle' 
spectrum with the same overall shape but with the acoustic 
oscillations erased. This is done mathematically by 



P(k) = Pun 



X + 



T nw (fc) x (1 



T lin (fc) 



(48) 



where Py m is linear matter power spectrum, T nw (fc) and 
Tn n (fc) are 'no-wiggle' and l i near t ransfer functions com- 
puted from lEisenstein fc Ht3 (|l998T ) and x = exp(— k 2 a 2 ) 
with a = 7 hT 1 Mpc chosen to fit the BAO suppression seen 
in their N-body simulations. 

The P(k) is then corrected for non-linear gravita- 
tional collapse using the HALOFIT software. The final 
P(k) is then transformed to £(r) using Eq. [19] Although 
the scale-dependent redshift-distortion an d halo bias cor- 
rectio n is weak at these scales, we follow lEisenstein et al.l 
(|2005h and multiply the correlation function by the square 



© 2011 RAS, MNRAS 000,m[26l 



22 U. Sawangwit et al. 



10 



r (h" 1 Mpc) 



100 



0.100 



0.001 r 



O w com (9) 
— - w com (9)-0.001 

w c (())-().()(! 13 

non-linear ACDM 

Projected Eisenstein ^(r) 

ACDM+low-1 excess 




10 



100 



8/ arcmin 



Figure 22. The combined w(6) (open circles) compared to the projections of non-l inear ACDM model pl us mode coupling, scale- 
dependent redshift-space distortion and halo bias correction (cyan solid line) and the lEisenstein et alj l|2005f ) £(s) (blue dashed line). 
The red dot-dashed line is the ACDM model plus low-Z power excess (" i|6.2.2t . The dash-dot-dottcd and dotted line shows the effect of 
sub tracting the data by 0.0 01 and 0.0015, respectively. The £(r) models used in the w(8) projec tion are given as an ins et together with 
the lEisenstein et al.l d2005Tl measurement (diamonds). Here, the same symbols are used for the lEisenstein et al. I J2005T ) and non-linear 
ACDM £(s) models as for the w(9) models above. 



of 1 + 0.06/ [1 + (0.06r) 6 ] (solid line in their Fig. 5), again 
chosen to fit what is seen in the N-body simulations. Such 
a correction is small at the BAO scale, only sub-percent at 
r > 25 ft -1 Mpc and increases to « 10% at 10 ft -1 Mpc. We 
then correct for the linear redshift-space distortion, the £(s) 
amplitude is enhanced re lative to the real-space correlation 
function, f(r), such that (|Kaiserlll987l ) 



2P 
3 



6' 2 



€(r) 



(49) 



Here, we assume /3 = 0.45 for these LRG samples (see 
iRoss et"al] l200?f l . The final f(s) model prediction with 
galaxy bias b = 2.09 for SDSS-L RG (see 34.211 is shown ( cyan 
solid line) in the inset of Fig. [22] lEisenstein et al.l (|2005h find 
this model to be a good fit to their £(s) data with the best- 
fit x 2 = 16.1 on 17 de grees of freedom for a particular set 
of cosmological parameters given above. We then computed 
w{9) from the £(r) derived above via Eg. 1141 Although the 
model (cyan solid line in Fig. 12211 was found to be consistent 
with the LRG £(s), it is inconsistent with our w(9) mea- 
surement, especially at r > 60 ft -1 Mpc or 8 > 120'. The 
uber-cal AAf2 w(9) between 40' — 400' (corresponding to 
20 < r < 200 /i" 1 Mpc) are incompatible with the model at 
99.8 percent level (x 2 =39.3 over 18-1 d.o.f with covariance 
matrix) . We note that this rejection may be associated with 
the apparent clustering excess at 6 > 200', which still could 
be subject to systematics. 

Next, we compare our w(9) to the best estimate of £(s ) 
at the BAO scale as measured by lEisenstein et al.l (|2005l ) . 



Although these measurements may have been superseded by 
DR7 spectrosopic LRG clustering analyses based on larger 
samples, these more recent estimat es are usually in reason- 
able agreement with the results of lEisenstein et all ll2005h. 



whether they are in c orrelation function (Martine z et al 



2009; iKazin et al.ll2010l ) or power spectrum (|Percival et al 



2010h form. For our comparison, we t hus simply make a poly - 



nomial fit to the best estimate £(s) of lEisenstein et all (|2005l ) 
(blue dashed line in the inset of Fig. I22p . The polynomial- fit 
£(s) is Kaiser de-boosted (Eq. I49|l to give £(r) by assuming 
P — 0.45. The £(r) is then corrected for the linear growth 
between z = 0.35 and z — 0.68 which reduces the amplitude 
by ~ 30%. The resulting model has similar amplitude with 
the expected AAf2-LRG £(r) because the SDSS and AAf2- 
LRG linear biases are coincidentally the same (see H4.2p . The 
model is then projected to w(9) using Eg. 1141 and is shown 
as a blue dashed line in Fig. 1221 Our result appears to be 
in good agreement with the model up to w 120 h~ l Mpc 
given statistical uncertainties in our measurement and the 
£(s) data. Beyond ~ 120 hT 1 Mpc, our w(8) shows a higher 
clustering amplitude as noted above. 

Summarising, the w CO m result appears consiste n t with 
the w(6) prediction based on the lEisenstein et al.l (|2005l ) 
best estimate of £(s) (at least out to w 120 Mpc) but not 
with the prediction based on the flat ACDM model due to 
the apparent large-scale clustering excess i n the w(9). This 
means that given the size of error bars of the lEisenstein et al.l 
(|2005h result, the ACDM model is quite compatible with 
the £(s) data but given the much smaller statistical error on 
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w{6), in this case our measurements are inconsistent with 
the ACDM model. While the feature observed at « 300' 
persists in most of the systematic tests we performed on the 
AAfi samples (i ]6.1[) . a few of these tests, e.g. dust extinc- 
tion, indicate there is still the possibility that systematic 
errors are affecting the w(6). Therefore, if we now assume 
that the excess signal at at 150 h -1 Mpc is an indication of 
a systematic and subtract 0.001 to 0.0015, the level of the 
excess amplitude at this point in iu com (see Fig. I22|) . we ob- 
tain the w{9) results as shown by the dash-dot-dotted and 
dotted lines. These two lines now bracket the flat ACDM 
result. Thus the issue of the disagreement between the w(9) 
result and the ACDM model seems to rest on the reality of 
the apparent clustering excess at large scales. 



Thomas et al. (2011) o 

best-fit ACDM 

AA£2-LRG (ACDM) - - - 
low-1 excess 




6.2.2 low-i power excess? 



Recently, IThomas et all (|2010|) ( 



also 



50 100 150 

Multipole moment (I) 



200 



Padmanabhan et alj 20071 ; Blake et al. 1 120071 ; IThomas" et al.l 

20 111 ) has also found a significant excess in their angular 



power spectrum, Ci, at the low multipoles relative to the 
best-fit ACDM models. They used photometric-redshift 
catalogues of the LRGs at z ~ 0.5 similar to our 2SLAQ 
sample. The most significant (« 4cr) low-^ power excess is 
observed in the highest redshift bin. The author carried out 
various systematic checks and found no indication of such 
an effect. While the clustering excess only affects Ct at 
multipoles smaller than the acoustic oscillations in Fourier 
space, unfortunately in configuration space the effect is 
expected on a wider range of scales and could affect our 
w(ff) BAO measurement. 

Fig. [231 shows the IThomas et all (120101 ) Ce and the ex- 
cess power model used to predict the AAf2 w( 8). For fur- 
ther d etail of the low-^ power excess model, see ISawangwitl 
|201ll ). The resulting w(6) model with the amplitude nor- 
malised to fit the data at = 40' — 400' is shown as the red 
dot-dashed line in Fig. 1221 The model appears to be consis- 
tent with our lOcom (13 percent confidence level, \ 2 = 23.6 
over 18 - 1 d.o.f) at r « 20 - 200 h' 1 Mpc. Therefore we 
note that the excess clustering signal in our w(8) is in good 
agreem ent with that observed in the Ci by IThomas et al.l 
|20ld ). The fact that the excess power in the Ce taken the 
form of an £ ~ 10 spike, suggests that this excess in w(8) 
is due to something other than acoustic oscillations in the 
power spectrum. 

We note that evidence for a large-scale (> 150 h^ 1 Mpc) 
correlation function excess has also been det e cted i n the 
NVSS ra d io so urce survey by iBlake fc WaH (|2002l ) and 
IXia et alj (|2010h . We have compared our results with theirs 
and find that our correlation function shows a similar shape 
but a factor of 2-3 x lower amplitude. If the excess clustering 
signal observ ed here is real th en it could be evidence for non- 
Gaussianity (|Xia et alj|2010l ) or for the gauge dep endence of 
the matter pow er spectrum on the largest scales (|Linll200ll ; 
IYoo et al]|2009l ) . But until this feature is detected in an inde- 
pendent galaxy dataset, there will always be the possibility 
that it is caused by some unknown systematics. Certainly, 
if the ACDM model were correct then we would have to 
conclude that this excess was caused by systematics at the 
level of Aw « 0.001 - 0.0015 in the photometric AAfi-LRG 
sample. 



Figure 23. The angular power spectr um of the 0.6 < z < 0.65 
McgaZ-DR7 LRG llThomas et al.l201ll) with significant power ex- 
cess at low multipoles (diamonds). The low-£ power excess plus 
the best-fit ACDM model of the 0.6 < z < 0.65 Ce (blue dot- 
dashed and solid lines) is extrapolated to the AAS1 redshift range 
(red dot-dashed and dashed lines). 



7 SUMMARY AND CONCLUSIONS 

We have presented here a new and detailed analysis of the 
angular correlation function of the Luminous Red Galaxies 
extracted from the SDSS DR5 photometric catalogue. All 
the necessary information for inferring the spatial cluster- 
ing is obtained and calibrated through redshift surveys of 
sample subsets. Our conclusions are as follows; 

(i) We measured the angular correlation function of the 
LRGs at three different redshifts, namely 0.35, 0.55 and 0.68 
and found the results to be well aproximated by power-laws 
at small and intermediate scales. 

(ii) With the large samples in terms of the numbers of 
objects and volume cover by the data, we see the deviation 
from the canonical single power-law at high significance. 

(iii) The data are better fitted by a double power-law 
where the large-scale ( > 1 — 2 h~ l Mpc) slope is equal to 
that of the conventional single power-law, i.e. 7 ~ 1.8. 

(iv) The form of the angular correlation functions at large 
scales are consistent with the expectation of the linear per- 
turbation theory in the flat standard ACDM Universe. 

(v) The LRG linear bias is high, b g ~ 2.0, as expected 
for massive luminous early-type galaxies and the clustering 
strength is found to be strongly linked to the sample intrinsic 
brightness. 

(vi) The best-fit HOD models suggest that these LRGs re- 
side in the masssive dark matter haloes, 10 13 — lO 14 ft _1 M0, 
and are typically central galaxies in their dark matter halo 
hosts, with the satellite fraction no more than 10 percent. 

(vii) The clustering evolution at intermediate scales (1 < 
r < 20 h~ l Mpc) is remarkably slow and may be approx- 
imately explained by a long-lived model or even a no- 
evolution model. The long-lived model may be in line with 
the observed passive evolution of the LRG luminosity func- 
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tion, consistent with a constant comoving LRG space den- 
sity in this redshift range. This latter conclusion would also 
apply in the case that the no-evolution (comoving) model 
were found to fit better but in this case the observations 
may require a signif icantly higher b ias. 

(viii) Using the lLacev fc Cole] (| 19931 ) framework, our 
Mdmh(z) measurements are well fitted by the model where 
halo mass is grown via merging of progenitors with masses 
of « 1.4 x 10 13 h,- 1 A/ Q and « 2.3 x 10 13 M from z = 1, 
for haloes that typically host L ^ 2L* and ^ 3L* galaxies, 
respectively. We found that these dark matter haloes have 
tripled their masses over the last half of cosmic time (al- 
though see the caveat given at the end of H5.1.2p whereas it 
has been claimed that t he LRG stellar m asses have grown 
by less than 50 percent (|Cool et alj|2008h . 

(ix) At small scales (r < 1 h^ 1 Mpc) the clustering evo- 
lution appears slightly faster at fixed luminosity and the 
clustering increases towards lower redshift, consistent with 
a virialised clustering model. Since our virialised model as- 
sumes a constant comoving LRG space density, a combina- 
tion of this stable clustering model at small scales and the 
long-lived model at intermediate scales could be consistent 
with the idea that merging of LRGs may not change the 
LRG space density significantly out to z ~ 0.7. 

(x) However, the evolution based on HOD and the ACDM 
halo merging framework requires that ~ 2 — 3 percent/Gyr 
of the LRGs merge with each other in order to explain the 
s mall-scale clustering evol ution, consis t ent w ith the results 
of I White et all (|2007T ) and I Wake et~al] (|200St ). 

(xi) In our AAS7-LRG result we find a BAO peak at a level 
consistent with the best estimate of £(s) obtained by Eisen- 
stein et al (2005). But, given the small size of our statistical 
errors, these results deviate significantly, ~ 4ct, from the 
standard ACDM prediction because of an apparent large- 
scale clustering excess. 

(xii) The excess clustering signal generally persists after 
a series of systematic tests we performed. However, a few 
of these tests did change the feature somewhat, suggesting 
that it could still be caused by some unknown systematic 
effects. 

(xiii) If the ACDM model were correct then we would 
have to conclude that this excess was caused by systematics 
at the level of Aw « 0.001-0.0015 in the photometric AAfi- 
LRG sample. 

(xiv) Otherwise, the excess signal in our w(0) relative to 
the standard ACDM model appears to be in good agreement 
with the Ct power excess at low I observed by other authors 
who used photo-z LRG samples atz« 0.5. 

(xv) If real, the large-scale clustering excess may be in- 
terpreted as an evidence for a non-standard cosmological 
model, e.g. primordial non-gaussianity or general relativis- 
tic effects. However, more, independent, data is required to 
check the reality of this clustering excess. 
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APPENDIX A: ANGULAR CORRELATION 
FUNCTIONS AND COVARIANCE MATRICES 

At the referee's request, we tabulate the angular correlation 
functions (Table |A"T|) measured from the three photometric 
LRG samples studied in this paper. The full covariance ma- 
trices in the form of correlation coefficients are also shown 
in Fig. |A"T1 



Table Al. The measured angular correlation functions for the 
SDSS, 2SLAQ and AAf2-LRG and their la JK errors. 
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Figure Al. The correlation coefficients, r.ij, out to very large angular separations. These are derived from the covariance matrices (Eq. 
I13H via 96 jackknife re-sampling fields. Three panels show r^j for SDSS, 2SLAQ and AAf2 -LRG samples from left to right. 
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